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ABSTRACT 

We have compared the 2001 XMM-Newton spectra of the stellar mass black hole binary 
XTE J1650-500 and the active galaxy MCG-6-30-15, focusing on the broad, excess emis- 
sion features at ^4-7 keV displayed by both sources. Such features are frequently observed 
in both low mass X-ray binaries and active galactic nuclei. For the former case it is gener- 
ally accepted that the excess arises due to iron emission, but there is some controversy over 
whether their width is partially enhanced by instrumental processes, and hence also over the 
intrinsic broadening mechanism. Meanwhile, in the latter case, the origin of this feature is 
still subject to debate; physically motivated reflection and absorption interpretations are both 
able to reproduce the observed spectra. In this work we make use of the contemporaneous 
BeppoSAX data to demonstrate that the breadth of the excess observed in XTE J 1650-500 is 
astrophysical rather than instrumental, and proceed to highlight the similarity of the excesses 
present in this source and MCG-6-30-15. Both optically thick accretion discs and optically 
thin coronae, which in combination naturally give rise to relativistically-broadened iron lines 
when the disc extends close to the black hole, are commonly observed in both class of object. 
The simplest solution is that the broad emission features present arise from a common pro- 
cess, which we argue must be reflection from the inner regions of an accretion disc around a 
rapidly rotating black hole; for XTE J1650-500 we find spin constraints of 0.84 ^ a* ^ 0.98 
at the 90 per cent confidence level. Other interpretations proposed for AGN add potentially 
unnecessary complexities to the theoretical framework of accretion in strong gravity. 

Key words: X-rays: binaries - Galaxies: active - Black hole physics - Individual Sources: 
XTE J1650-500, MCG-6-30-15 



1 INTRODUCTION 

Simple models for the accretion of material onto a black hole 
(of any mass) usually include two physically distinct components. 
The first is an optically th ick, geometrically thin accretion disc 
jShakura & Sunvaevlll973l) formed by material which is gravita- 
tionally bound to the black hole, but has sufficient angular momen- 
tum to av oid falling dire ctly onto the black hole in the manner dis- 
cussed by iBondJ l ll952h . The second is a corona of hot electrons, 
analogous to that seen around the Sun (although the exact geome- 
try and energy distribution of these electrons is not yet well deter- 
mined, see e.e. ICoDpil 19991) . populated by both ther mal and proba- 
bly non-thermal electrons jMalzac & Belmont! 20091) . The emission 
processes relevant to each of these components are relatively well 
understood. Viscous interactions within the disc heat the material 
and thermalise the emission, allowing the disc to be treated as a se- 
ries of black bodies of increasing temperature with decreasing ra- 
dius. Some fraction of the thermal emission from the disc is then re- 
processed by the coronal electrons, and Compton up-scattered into 
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a powerlaw-like emission component jHaardt & Maraschilll99lh . 
with a high energy cut-off determined by the temperature of the 
thermal electrons, and a low energy roll-over determined by the 
temperature of the accretion disc. 

The temperature at which the thermal emission from the disc 
is expected to be observed depends on the mass of the black hole. 
For the black holes associated with X-ray Binaries (XRBs; Mbh 
~ 10 Mq) the discs are relatively hot and the thermal emission 
should be observed in the X-ray band at ~0.2-1.0keV, while the 
discs around the much larger black holes associated with active 
galactic nuclei (AGNs; A/bh > lO'^ Mq) are much cooler and 
should be observed in the ultraviolet (UV). However, in both cases 
the temperature of the coronal electrons is so high (~100keV) that 
the Comptonised emission is always observed in X-rays. Observa- 
tions of both stellar and supermassive black holes in the UV and 
X-ray bands have confirmed the presenc e of these emission com- 
ponen ts; fo r XRBs see recent reviews bv lRemillard & McClintocfl 
booeh and lDone et all ( l2007h . for AGNs see lElvisI feOld) and ref- 
erences therein. 

With the onset of high-resolution CCD spectroscopy it has 
become clear that the broadband X-ray spectra of AGN are in 
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many cases far more complex than the simple powerlaw contin- 
uum expected to be produced by the black hole corona. Addi- 
tional features arise which can, depending on the intrinsic contin- 
uum adopted, be interpreted as originating from either emission 
or absorption processes, and viable, physically motivated models 
have been proposed for each case. In emission interpretations the 
two most prominent features are a broad and smooth excess below 
^l.OkeV, commonly referred to as the 'soft excess', and a sec- 
ond broad excess over ~4.0-8.0 keV, while in absorption models 
these are instead interpreted as a broad absorption trough spanning 
~1.0-5.0keV. The energy of these features, however interpreted, 
remains remarkably constant, despite the observed AGN host- 
ing black holes spanning many orders of magnitud e in mass (see 
e.g. iGierlihski & Don^ |2004 ICrummv et al.l |200€ iNandra et al.l 



l2007l iMiniutti et al.ll2009l . iBrenneman & Revnoldsl l2009h . which 
strongly suggests that if they are due to additional emission they 
must also originate through atomic processes (absorption being an 
atomic process by nature). 

One of the successful emission models is the disc reflec- 
tion/light bending interpretation which assumes that the corona 
is compact (or centrally concentrated) and emits some fraction 
of its flux towards the accretion disc. This is reprocessed into 
an additional 'reflected' emission component via backscattering, 
fluorescence, an d Compton scattering in teractions within the sur- 
face of the disc dOeorge & Fabianlll99ll) . This reflection compo- 
nent includes discrete atomic features imprinted by the disc, the 
most prominent of which is usually the neutral iron Kq emis- 
sion line at 6.4 keV. However, the strong gravity close to a black 
hole broadens and skews thes e features via Doppler shifts, b eam- 
ing and gravitational redshifts l lFabian et alJl989l : lLaodll99ll) . The 
broad excess over ~4.0-8.0keV is thus attributed to the iron Ka 
emission line broadened by such relativistic effects, and the large 
number of emission lines from iron and lower-mass elements in 
the ~0.5-1.0keV range are blended together into a broad and 
smooth emission component, similar to that observed. The strong 
gravity also focuses some additional amount of the hard emission 
down onto the disc through gravitational light bending. This phe- 
nomenon causes the spectral variability often observed; if the re- 
gion in which the Comptonised continuum is produced varies in 
distance from the central black hole, the fraction of the coronal flux 
foc ussed onto the disc, and hence that observed, also varies (see 
e.g. IMiniu tti & Fab iarj|2004h. For a recent review on X-ray reflec- 
tion, seelfabian & Rosj 1 I2OIOI) . 

An alternative interpretation is that the observed deviations 
from the powerlaw continuum arise as a result of complex absorp- 
tion due to structures of intervening clouds and infalling/outflowing 
material. Combinations of these structures can lead to complex 
combinations of neutral, partially and fully ionised material. In- 
deed, in a number of AGN there are clear, undeniable signatures of 
absorption due to partially ionised material (see e.g. iBlustin et al.l 
I2OO5I) ; these 'warm' absorbers are usually associated with mildly 
outflowing material from the dusty torus or clouds in the broad 
line region. In addition, there are also detections of absorption by 
highly ionised material with mildly relativistic blueshifts claimed, 
which are interpreted as evidence for outflowing disc winds (see 
iTombesi etid]|20I(l iReeves et alj|2009l. IPounds & ReevesI l2009l 
and references therein), although it ma y be necessary to treat a 
number of these claims with caution (see lVaughan & Uttlevl200i) . 
However, in order to reconcile the large extent of the apparent ab- 
sorption trough in energy (up to ~5 keV) with its relatively shallow 
appearance, it is usually necessary to invoke absorption by partially 
ionised material that only covers some fraction of the X-ray source. 



With appropriate combinations of these various types of absorp- 
tion the observed spectra can be reproduced, and spectral variability 
may be explained via changes in the covering fraction, and/or col- 
umn density and ionisation of the absorbing mate rial. For a review 
of the absorption processes relevant to AGN see [Turner & Milled 
( l2009l) . 

Even with the quality and quantity of data provided by the 
latest generation of X-ray observatory, such as XMM-Newton and 
Chandra, it has proven very difficult to statistically distinguish be- 
tween these interpretations for AGN, a c ase in point being the 
highly variable A GN MCG-6-30-15 (see iFabian etal]|2002l and 
[Miller et al.[|2008| , who construct statistically acceptable reflection 
and absorption dominated models respectively). In principle, the si- 
multaneous coverage of the ~0.5-50.0 keV energy range possible 
with the Suzaku satellite may be able to distinguish between disc 
reflection and pure absorption interpretations in some cases. Re- 
flection models predict the presence of an additional, broad emis- 
sion feature at ~20-30 keV, referred to as the Compton Hump as 
it arises due to the interplay of Compton up-scattering of low en- 
ergy photons and photoelectric absorption of high energy photons 
within the reflecting medium. Features consistent wi th the Comp- 
ton h ump are fairly commonly seen in AGN (see e.g. lWalton et al.l 
I2OIOI) . However, owing to the relatively poor quality data currently 
available above lOkeV in many cases, these high energy features 
can also be consistent with the presence of a partially covering, 
Compton thick absorber (see e.g. iTumer et al.ll2009l : iReeves et al.l 
\2QQ% . 

One of the strongest individual cases in favour of the reflection 
interpretation is the narrow-line Seyfert I galaxy I H 0707-495, in 
which reverberation betwe en the powerlaw and reflection compo - 
nents has been detec ted (see lFabian et alj2009llZoghbi et al.l201(lh . 
iMiller et al.l l l20I0bh argue such apparent reverberation can still be 
reproduced by distant reflection from the winds required by an ab- 
sorption dominated interpretation. However, no absorption model 
that simultaneously reproduces the reverberation and the observed 
energy spectrum is offered. In addition, it is worth noting that soft 
excesses are observed in AGN which display no obvious indication 
of any absorption over that exp ected due to the Galactic interstel- 
lar m edium (ISM), e. g. Ark 120 jVaughan e t al.l2004lNardim et all 
I2OIII) and Ton SI 80 l lVaughan et al. 2002). 

Interestingly, high quality X-ray spectroscopy of XRBs has 
also revealed the presence of broad excesses over the Comptonised 
continuum in the ~4.0-8.0keV energy range, which are again 
frequently interpreted as r elativistically broadened iron emissio n 
lines, e.g. XTE J1650-500 dMiller et alj2002l;l> 
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GX 339-4 l lReis et al.ll2008l). GRS 1 9I5-HI05 ( iBlum et al.ll2009l). 
Swift J1753.5-0127 ( iReis et alj|2009l) . XTE J1752-223 l lReis et all 
I2OIOI) . In addition, high energy features consistent with the Comp- 
ton H ump are also frequently observed in BHBs (see e.g. lReis et al.l 
bold) , supporting the theory that X-ray reflection is present in 
these sources. Studying XRBs has the advantage that they are of- 
ten brighter in flux than their AGN counterparts, so the presence 
of absorption due to outflowing winds, etc. can be tested in greater 
detail, and their contribution, if any, to the observed spectrum de- 
termined with much greater confidence. 

However, the presence of relativistically broadened iron lines 
in XRBs is not without its own controversy. Although it is generally 
accepted that these features do arise due to iron K a emission, and 
that this emission has some finite width, there is some debate over 
how broad the lines actually are, in particular whether instrumental 
effects act to artificially enhance the apparent width and skew of the 
line profile, and therefore over the mechanism that broadens them. 
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This is entwined with the debate over whether the accretion disc 
extends in to the inner-most stable circular orbit (ISCO). While it 
is generally accepted that this is the case at high accretion rates, and 
that at some low accretion rate the optically thick disc most likely 
begins to truncate at a larger radius within which a hot, advection 
dominated accretion flow arises, the accretion rate at which trun- 
cation occurs is still under dispute. Some authors, e.g. iDone et al] 
( |2007|) . argue that the disc begins to truncate during the transition 
bet ween the h i gh/sof t and low/hard accretion states, while others, 
e.g. lReis et al 1( l2010l) . argue that the disc remains at the ISCO until 
much lower accretion rates (f or a review of the 'standard' accr etion 
states displayed by BHBs see lRemillard & McClintockll2006l) . 

If, as is likely, the line is produced in the disc, then as the 
disc recedes the line profile should evolve and become narrower. 
However, BHBs in the low/hard state still appear to di splay sim- 
ilarly broad excesses to other accretion states (see e.g. iReis et al] 
I2OIOI) . iLiu erai] ( [201 l|) proposed a modification to the truncated 
disc model, in which the outer disc truncates, but some inner por- 
tion remains, resulting in a disc with a hole in it. Such a geome- 
try would allow the truncated disc model and the broad iron line 
observations to be reconciled without invoking a separate medium 
that produces the line. However, through analysis of the variabil- 
ity displayed by both the therma l and Comptonised compon ents of 
GX 339-4 in the low/hard state. Iwilkinson & UttlevI | |2009|) show 
that in order for the variability to be driven by fluctuations in the 
mass accretion rate, as is generally thought to be the case, the disc 
must extend uninterrupted to at least 20i?G, placing strict con- 
straints on the location and size of any hole which render this sce- 
nario unlikely. The truncated disc interpretation therefore requires 
that the dominant line broadening mechanism is some process other 
than Doppler/gravitational effects, as this interpretation frequently 
requires that the disc extends close to the black hole. 

An alternative mechanism that has been invoked to explain 
broad line profiles is Compton broadening, in which Compton scat- 
tering of line photons modifies their energy (see e . g. .Czernv et al 
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20091) . This process broadens and, depending on the conditions of 
the Comptonising region, can even skew the line profile. Some con- 
tribution from Compton broadening is unavoidable for any realistic 
method of line production, as the line emission will necessarily be 
Compton scattered as it escapes from the medium in which it is 
produced. If the line is associated with the disc, Compton broad- 
ening is expected to be more significant process for BHBs than for 
their AGN counterpar ts owing to their relative disc temperatures 
( lRoss&Fab"iaij|2007h . The challenge then becomes determining 
the relative contribution to the line profile of this process and the 
relativistic effects discussed above. However, it has also been pro- 
posed that the line does not originate in the disc, but is instead pro- 
duced via recombination in the corona or in a disc wind. These 
scenarios require Compton broadening to play a prominent role in 
determining the line profile. We note that, given the changes in the 
corona and the changes in mass inflow and outflow rate typically 
observed during the evolution of BHB outbursts, evolution in the 
line profile would also be expected for these latter scenarios. 

Here, we present an analysis of the high quality X-ray spec- 
tra of the XRB XTE J1650-500, which displays a broad line pro- 
file. Both relativistic and Compton broadening are examined, and 
we demonstrate that the former scenario is strongly preferred. We 
then proceed to compare the high energy spectra of XTE J1650-500 
and the AGN MCG-6-30-15, and demonstrate the similarity of the 
features identified as broad iron lines in each case. This work is 
structured as follows: section |2] gives details on the data reduction. 



section [3] presents our spectral analysis, highlighting the similarity 
of the two line profiles, and finally section |4] presents our discus- 
sion, in which we use this comparison as the basis for a simple, 
logical argument in favour of the reflection interpretation for AGN. 



2 OBSERVATIONS AND DATA REDUCTION 

XTE J 1650-5 00 and MCG-6-30 -15 were both observed with 
XMM-Newton jjansen et al]|200ll) in 2001 September and August 
respectively. In addition, XTE J 165 0-500 was observed three times 
by BeppoSAX l lBoella et al.lll997ij) during the same outburst. The 
first BeppoSAX and the XMM-Newton observations of XTE JI650- 
500, separated by about a day, caught the source in a rising in- 
termediate state, i.e. during the transition between the low/hard 
and high/soft states in which the thermal and Comptonised emis- 
sion components are both strong. The latter two BeppoSAX ob- 
servations, taken further along the evolution of the outburst, found 
the source in a more traditional high/soft state. These observations 
were selected because they represent some of the best quality X- 
ray data available for XRBs and AGNs, as the data obtained have 
good photon statistics, the XMM-Newton observation provides low 
energy coverage and the BeppoSAX data provides a view of the 
emission from XTE J 1650-500 at ~6keV with an independent de- 
tector. Here we detail the data reduction procedure adopted for each 
source. In both cases the XMM-Newton reduction was carried out 
with the XMM-Newton Science Analysis System (SAS vlO.0.0). 



2.1 XTE J1650-500 

The 2001 observation of XT E J 165 0-5 00 wa s taken on September 
13th with the EPIC-pn ( iStruder et al.ii200li) CCD in burst mode, 
so despite its extremely high flux, the spectrum obtained with this 
detector is not expected to suffer from pile-ujj^. The data reduction 
procedure adopted for this source follows closely that outlined in 
iKirsch et al.l 1 120061) . The EPIC-pn observation data files were pro- 
cessed using EPCHAIN to produce calibrated event lists from which 
spectral and timing products may be extracted. Burst mode is simi- 
lar to Timing mode and operates with extremely fast readout times; 
200 lines of CCD pixels are fast-shifted in 14.4 ps as the source 
data is being recorded. This readout process leads to loss of spa- 
tial information in the shift (RAWY/CCD line number) direction, 
such information is only available in the direction on the CCD per- 
pendicular to the readout (RAWX/CCD column number). Despite 
the ex tremely fast readout employed in burst mode, Kirsch et al] 
( I2OO6I) found that events with RAWY > 140 contain piled-up data, 
so we extracted the source spectrum from the region bounded by 
s; RAWY ^ 140 and 27 ^ RAWX < 46. In addition, we also 
only included single and double events, and excluded border pixel 
events. As XTE J1650-500 was so bright during this observation, 
there was no region on the operational CCD free of source counts, 
so we did not extract a background spectrum following the rec- 
ommendation of the XMM-Newton burst mode reduction guid^fl 
Given the high source countrate (see below), the contribution of 
background counts should be negligible. The redistribution matri- 
ces and auxiliary response files were generated with RMFGEN and 



^ Pile-up refers to the situation in which more than one photon is inci- 
dent on a CCD pixel within a single readout time. These are incorrectly 
registered as a single photo n at an a rtificially high energy. For a detailed 
description of this effect see lMillei~et al. 1 20 10a) and references therein. 
^ http://xmm.esac.esa.int/sas/current/documentation/threads/ 
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ARFGEN respectively. The spectrum obtained had a good exposure 
time of 685 s, and an average total countrate of 2504 ± 2 ct s~^. 
Finally, the spectrum was grouped using GRPPHA such that each 
energy bin contained a minimum of 25 counts, so that the probabil- 
ity distribution of counts within each bin can be considered Gaus- 
sian, and hence the use of the statistic is appropriate when per- 
forming spectral fits. The EPIC-MOS detectors were not operated 
for this observation, and we do not consider the spectrum obtained 
with the RGS instrument in this work, as it suffers quite severely 
from pile-up effects. 

Throughout this work, we also make use of the three Bep- 
poSAX observations taken during the course of the same outburst. 
These were obtained on 2001 Septembe r 11th, 21st and Oc tober 
3rd respectively. We reduced the MECS l lBoella etalj|l997bl) data 
for each observation following the prescription outlined in the Bep- 
poSAX cookbool0. In each case we defined a circular source region 
of radius 8', and a circular background region free of other sources 
of the same size, and use the appropriate September 1997 auxil- 
iary response and redistribution files. The good exposure times and 
average countrates obtained for each observation are 47, 68 and 
24 ks and 34.18 ±0.03, 42.18 ±0.03 and 32.57 ±0.03 ct s"^ re- 
spectively. Note that the MECS detector has a lower energy resolu- 
tion than EPIC-pn (~750 eV in comparison to ~ 150 eV at 6 keV). 



2.2 MCG-6-30-15 

The long 2001 observation of MCG-6-30-15 was taken in small 
window mode, and spanned three orbits. Data reduction was car- 
ried out separately for each orbit, largely according to the stan- 
dard prescription provided in the online guid^fl- In this work we 
focus on the EPIC-pn data for direct comparison with XTE J1650- 
500. The observation data files were processed using EPCHAIN to 
produce calibrated event lists. Spectra were produced for the 0.5- 
10.0 keV energy range selecting only single and double events us- 
ing XMMSELECT, and periods of h igh background were t reated 
according to the method outlined bv IPiconcelli et al] ( |2004|) . with 
the signal-to-noise ratio maximised for the 6.0-10.0 keV energy 
band, where the neutral and ionised iron transitions are expected 
to be observed. A circular source region of radius 37" was de- 
fined in order to include as many source counts as possible without 
including CCD chip edges, and a larger circular background re- 
gion of radius 54" was chosen in an area of the same CCD free of 
other sources. The redistribution matrices and auxiliary response 
files were generated with RMFGEN and ARFGEN. After perform- 
ing the data reduction separately for each of the orbits, the spec- 
tra were combined using the FTOOI0 ADDSPEC, which combines 
spectra in a response weighted manner to account for any dif- 
ferences there may be in e.g. the exposure times, etc. ADDSPEC 
also automatically combines the instrumental responses and back- 
ground spectra associated with the source spectra. The resulting 
spectrum had a good exposure time of 221 ks, and an average total 
count rate of 28.15 ±0.01 ct s~^. Finally, for the same reason as 
for XTE J1650-500, the spectrum was grouped using GRPPHA to 
have a minimum of 25 counts in each energy bin. 



^ http://heasarc.gsfc.nasa.gov/docs/sax/abc/saxabc/saxabc.html 
^ http://xmm.esac.esa.int/ 

^ http://heasarc.nasa.gov/ftools/ftools_menu.html 



3 SPECTRAL ANALYSIS 

Here we detail our spectral analysis of XTE J1650-500 and MCG- 
6-30-15. The approach for XTE J1650-500 is simple: we begin by 
modelling the continuum to highlight the broad line profile, be- 
fore modelling the line in a phenomenological fashion and finally 
moving on to a more physical treatment of the disc reflection. In 
addition, we also consider Compton broadening in some detail and 
find that strong relativistic effects are still required to explain the 
breadth of the line. We then construct a simple phenomenological 
model for the spectrum of MCG-6-30-15 in order to again high- 
light the line profile in this source, and its similarity to that of 
XTE J 1650- 500 (a full physical inte rpretation of MCG-6-30-15 is 
presented in lChiang & Fabianl201 ll . utilising all currently available 
data, and so will not be undertaken here). Spectral analysis is per- 
formed with XSPEC vl2.6.0f (Arnaud 1996), and quoted uncertain- 
ties on model parameters are the 90 per cent confidence limits for 
one parameter of interest determined by variation, unless stated 
otherwise. In both cases, due to calibration uncertainties possibly 
related to a deficient charge transfer inefficiency (CTI) correction at 
high count rates (see section U. 1.7l for a description of the effects of 
CTI), significant negative residuals were continually observed as- 
sociated with the instrumental silicon K and gold M edges at ~ 1.8 
and 2.2 keV respectively, hence the energy range 1.7-2.5 keV has 
been excluded from our anal ysis throughout thi s work. We adopt 
the solar abundances given in lWilms et al. 

3.1 XTEJ1650-500 

3.1.1 Continuum Modelling 

We model the continuum of XTE JI650-500 with the standard 
black hole binary (BHB) components, namely a multi-colour black 
body accretion disc (modelled with DISKBBi lMitsuda et al1ll984h 
which accounts for the thermal emission from the disc and domi- 
nates at soft energies, and a powerlaw, associated with Compton- 
isation of the thermal emission by electrons in a corona, which 
dominates at hard energies. Both of these are modified by Galac- 
tic absorption, acc ounted for with the TBABS absorption code 
jWilms et al.ll2000l) . The column density of the ISM, as well as its 
Oxygen, Neon and iron abundances, were allowed to vary; all other 
elements were assumed to have solar abundances. The 4.0-8.0 keV 
data were not considered at this stage of the analysis, in order to 
minimize any contribution from the iron emission region to our de- 
termination of the continuum. This model was fit to the remaining 
data, and although statistically it is not formally acceptable, with 
xt = 1184/865, it provides a good representation of the general 
shape of the spectrum. The relatively poor fit statistic is mostly due 
to residual features at ~0.6-0.8 keV (see Fig.lTJ, which will be dis- 
cussed in section [3.1.51 the data closer to the 4.0-8.0 keV range, 
where the iron line is expected to be observed, is well modelled. 

3.1.2 Line Profil e 

We now include the 4.0-8.0 keV data in our analysis. Fig.[T](panel 
A) shows the data/model ratio of the full 0.5-10.0 keV spectrum for 
XTE J1650-500 to the previously determined continuum model. 
It is clear that there is a large, broad excess in the ratio spec- 
trum in this energy range, which we interpret as arising due to 
iron emission. Initially, we modelled the line profile with a Gaus- 
sian emission line component, with the rest-frame energy of the 
line constrained to be between the neutral iron Ka transition at 
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Table 1. Key continuum and iron line parameters obtained by phenomenologically modelling the broad line profiles of XTE J1650-500 and MCG-6-30-15 
with LAOR line models. Parameters marked with * have not been allowed to vary, and those with 'p' quoted as one of the limits have reached the relevant 
constraint imposed upon it. The rest frame energies of the iron emission lines were required to be in the range 6.40-6.97 keV. The temperatures quoted are of 
the DISKBB component for XTE J1650-500, and of the black body included to account for the soft excess in MCG-6-30-15. 
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6.4 keV and the iron Lya transition at 6.97 keV, as the d iscs in 
XRBs are likely to be highly ionised i IRoss & Fabiarill2007h . With 
the inclusion of the Gaussian line the model provides a fairly good 
fit to the data, with xt = 1936/1663. The soft residuals high- 
lighted previously are still present, but the line profile is well mod- 
elled. We find that the line energy is constrained to the lower end 
of the allowed energy range in contrast with the expectation for 
an ionised disc, with iScauss < 6.5 keV, and the profile is very 
broad, wither = l.nt°ol\ 

keV. The equivalent width obtained is 

EW = 350l}?° eV. 

If the iron line is primarily broadened by gravitational and 
Doppler effects, then its width implies the disc must extend to the 
regions of strong gravity close to the blac k hole. We r eplaced the 
Gaussian emission line with a LAOR2 line ( lLaoj|l99l]) , which ap- 
plies the relevant relativistic effects to a Gaussian emission line 
profile. The parameters of this component are the energy of the 
emission line, the inner and outer radii of the accretion disc in 
gravitational radii (Rg= GM/c^), the inclination of the disc, i, 
the indices of the radial emissivity profile, assumed to be a broken 
powerlaw of form: 



e(r) cx 



r > Rh 



(1) 



and the break radius of the emissivity profile, also in gravitational 
radii. Recent simulations have suggested the radial emissivity pro- 
file may indeed be more complex than a s i ngle powerlaw, showing a 
break at ~6 Rg dWilkins & Fabiaij201lljMiniutti & Fabiaij2004h . 
The outer radius of the disc was taken as the maximum allowed 
by the model, 400 Rg (this parameter is rarely well constrained as 
even in the case of a Schwarzschild black hole the majority of the 
emission arises from the innermost regions of the disc). The rest 
of the line parameters were initially free to vary, although in the 
course of the modelling we found that i?br and Qoutcr could not 
be constrained. These were respectively set to i?br = GRg, after 
which the metric tends asymptotically towards the Schwarzschild 
case even for a maximally rotating Kerr black hole, and qouta = 3, 
as expected for a lamp-post illumination model. With the inclusion 
of the LAOR2 line, the model again provides a fairly good fit to the 
data, with xt = 1936/1661. The parameters and equivalent width 
(EW) obtained for the LAOR2 profile are given in Table[T]and the 
data/model ratio plot is shown in Fig. [T] (panel B). The energy of 
6.97 keV obtained for the iron line with this profile is more consis- 
tent with the expectation of a highly ionised disc. In Fig. [T] (panel 
C) we also show the ratio plot for the model including the LA0R2 
line, but with the line normalisation set to zero, to demonstrate that 
this component models the hard residuals well without significantly 
altering the continuum model determined previously. 
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Figure 1. Data/model ratio plots for XTE J1650-500. Panel A: ratio plot 
to the continuum model fit to the 0.5^.0 and 8.0-10.0 keV data (see sec- 
tion [TTTT); the broad iron line profile is clearly shown. Panel B: the ratio 
plot when a LAOR2 line is included in the model; the iron line is clearly 
well modelled. Panel C: ratio plot to the model including the LAOR2 line 
with the hne normalisation set to zero, to demonstrate that this component 
accounts for the hard residuals without drastically altering the continuum. 
Soft residual features over the ~0.6-0.8 keV energy range are seen through- 
out, these are discussed in section [3. 1.41 and section [3. 1.61 Panel D: ratio 
plot of the 2002 burst mode spectrum of the Crab Nebula modelled with 
an absorbed powerlaw (see section [3. 1.4) . Similar soft residuals are seen, 
implying these features possibly arise due to systematic uncertainties. The 
spectra have been rebinned for clarity. 



3.1.3 High Energy Comptonisation 

Although the high energy Comptonised continuum is often 
powerlaw-like in the XMM-Newton bandpass, we also investigated 
whether this approximation has any effect on the line profile. To 
test this, we also modelled the high energy co ntinuum with the 
more physical Comptonisation codes COMPTT ( lTitarchuk|[T994h 
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Table 2. Key iron line parameters obtained when modelling Comptonised 
continuum with the COMPTT and COMPPS codes. The line profiles have 
again been modelled with a LAOR2 component, in the same manner as 
with the simple powerlaw continuum. The values obtained show excellent 
agreement with those presented in Table[T). 



Model 


E 




•Jinner 


i 




(keV) 


(ifc) 






(deg) 


COMPTT 


-0.32 


1.5 ±0.1 


> 


6.6 




COMPPS 


_0.44 


1.5 ±0.1 


> 


8.7 


69^7 



and COMPPS jPoutanen & Svensson|[l996h . In both cases the tem- 
perature of the seed photon spectmm was set to that of the DISKBB 
component, and we assumed a standard slab geometry. These mod- 
els were applied in the same manner detailed previously {i.e. ini- 
tially excluding the 4.0-8.0 keV data), and the line profiles ob- 
tained examined. As shown in Fig. |2l use of these more sophis- 
ticated Comptonisation models does not lead to any modification 
of the line profile. With the addition of a LA0R2 line we obtain 
fits of equivalent quality to the case using the simpler powerlaw 
continuum in both cases 

{xl = 1933/1660 and 1941/1660 respec- 
tively). The parameters obtained for the line profile with each of the 
models are given in Table |2l and show excellent agreement within 
their statistical uncertainties. Therefore, we conclude that a simple 
powerlaw component provides an adequate representation of the 
Comptonised continuum in the bandpass considered. We will pro- 
ceed with the use of this approximation throughout the rest of this 
work. In addition, it is clear that any systematic uncertainties in- 
troduced through our choice of continuum model are negligible in 
comparison to the statistical parameter uncertainties obtained. 

3.1.4 Disc Reflection 

We now embark upon a more physical consideration of the rela- 
tivistic disc reflection scenario for XTE J 1650-500, and replace 
the DISKBB and LA OR components with the REFBHB model 
( IRoss & Fabian ''2007) convolved with the KDBLUR kernel, which 
applies the same relativistic effects included in the LAOR line model 
to any general input spectrum, but assumes a simpler single power- 
law emissivity profile of form e(r) oc r~'. REFBHB is a self consis- 
tent reflection model which intrinsically calculates the iron K-shell 
absorption and emission in the reflecting medium, correctly includ- 
ing effects such as Compton broadening, which, as emphasised pre- 
viously, is an i mportant process in BH Bs due to the high tempera- 
ture of the disc jRoss & Fabianll2007l) . It is calculated specifically 
for use with spectra from BHBs as it includes the thermal emission 
from the accretion disc. Its key parameters are the hydrogen num- 
ber density and temperature of the accretion disc (nu and kT), the 
photon index of the illuminating powerlaw continuum, the relative 
strength of the illumination with respect to the thermal emission at 
the surface of the disc {Fu/Fth), and the redshift of the reflecting 
medium. We assume the redshift is zero, while the photon index 
of the illuminating emission was required to be the same as that of 
the powerlaw continuum component. This model provides a fairly 
good representation of the data, with = 1832/1663, and the iron 
line profile is well modelled. The key parameters obtained are given 
in Table[3]as Model 1, and the relative contributions of the compo- 
nents are shown in Fig. [3] As suggested by the results obtained 
with the LA0R2 line profile in section [3.1.2| we again find that the 
inner radius of the disc is small, with Rin ~ 2Rg, although the 
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Figure 2. Data/model ratio plots for the XMM-Newton observation of 
XTE J 1650-500, focussing on the iron emission, after making use of vari- 
ous models for the high energy Comptonised continuum. The models inves- 
tigated are a simple powerlaw (top panel), and the COMPTT and COMPPS 
Comptonisation codes (middle and bottom panels respectively). The choice 
of model does not effect the line profile obtained. 



radius obtained here is slightly larger than with the simple LAOR2 
profile, owing to the self consistent inclusion of Compton broaden- 
ing in REFBHB. In addition, the value obtained for Fu/Fth suggests 
that the Comptonised emission contributes significantly more to the 
total flux than the thermal emission from the disc, which is not un- 
common for a source during the transition from the low/hard to the 
high/soft state jRemillard & McClintoc32006l) . 

In order to assess whether our choice of reflection model has 
any major influence over the key parameter values obtained, we 
replaced REFBHB with REFLIONX iRoss & Fabianji2ol)5l) . another 
self consistent reflection code. In particular, we are interested in the 
inner radius of the disc, which in turn provides information on the 
spin of the black hole and is key in determining the width of the 
line. The key parameters of REFLIONX are the iron abundance and 
ionisation parameter (defined as ^ = L/nR^, where L is the ion- 
ising luminosity of the source, n is the density of the absorber and 
R is its distance from the source) of the reflecting medium, and the 
photon index of the ionising continuum. We again required the pho- 
ton index parameters of both the reflection model and the powerlaw 
continuum to be the same, and the iron abundance was set to the so- 
lar value. As expected, the reflecting medium is found to be highly 
ionised, with log^ = 2.8. However, although REFLIONX includes 
treatment of all the same processes as REFBHB, it is calculated for 
use with X-ray spectra from active galaxies, and does not include 
the direct thermal emission from the disc in the correct fashion for 
BHBs. In order to account for this emission we again included a 
DISKBB component. This alternative reflection model is also able 
to reproduce the data fairly well, with = 1848/1662. The inner 
radius obtained for the disc here is 7?in = LSllgJ^, even smaller 
than in the REFBHB case. This decrease can be understood in terms 
of the differing situations for which the two reflection models are 
calculated. We stress again that REFLIONX is calculated for use 
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Figure 3. The basic disc reflection model (Model 1) for XTE J 1650-500 
presented in i]3.1.4l and the relative contributions of the components in- 
cluded. The solid black line shows the total model, while the blue shows 
the powerlaw continuum and the red shows the blurred reflection compo- 
nent (which includes the thermal contribution from the disc). The spectra 
have been rebinned for clarity. 

with active galaxies, and therefore assumes a much cooler temper- 
ature and lower density for the accretion disc. This in turn reduces 
the contribution from Compton broadening, hence additional grav- 
itational broadening is required to reproduce the line profile, and 
a smaller inner radius is obtained. For this reason, the inner radius 
obtained with REFBHB is more appropriate as it includes the correct 
amount of Compton broadening for XRBs. 

It is worth briefly noting that the neutral column density 
for the ISM obtained is consistently larger tha n that presented in 
some of the previ ous work on this source {e.g. iMiller et alJllOOi. 
iMiller et"al]|2009l) . possibly owing to the inclusion of data in the 
range 0.5-0.7 keV in our analysis. If we restrict ourselves to a 
consideration of th e same 0.7-10.0 keV energy range analysed in 
IMiller et"al] | |2009|) . we obtain equivalently acceptable fits with a 
column density consistent with theirs. However, in doing so the 
abundances of Oxygen, Neon and iron obtained are all ~2 times 
solar. The difference in column obtained is therefore most likely 
to be due to the updated absorption models, solar abundances and 
cross-sections used in this work. In addition, none of the models 
presented required the energy of the neutral Oxygen recombina- 
tion edg e to be shifted away from the expected value, as originally 
found in lMiller et all j2002f) . 

3. 1.5 Soft Residuals 

In both cases the residuals at ~0.6-0.8 keV are still present. In or- 
der to test whether these features are real, or arise due to prob- 
lems with the burst mode calibration at low energies, we also 
briefly analyse the 2002 XMM-Newton burst mode observations 
of the Crab Nebula, obtained approximately six months after the 
XTE 11650-500 data presented. Data reduction was performed sep- 
arately for each of the 3 orbits following the prescription outlined 
in section 12.11 and the resulting spectra were combined again in 
a response weighted manner using ADDSPEC. The final spectrum 
obtained was modelled with an absorbed powerlaw, again using 
TBABS to model the Galactic absorption. We adopt t he abundances 
for the ISM in the direction of the crab reporte d by| Kaastra et al.l 
( l2009h (here adopting the solar abundances of iLodd ers 200^ as 
used in that work). The best fit column density and photon index 



are Nh^ 4.03 ± 0.02 x 10^^ atom cm"^ and T = 2.070 ± 0.004, 
although the fit qu ality is relatively poor, with xi = 2362/1743. 
iKirsch et al.l ilOOdb also obtained a fairly poor quality fit with a 
similar photon index to that obtained here. The data/model ratio 
plot obtained with this model is shown in Fig. [T] (Panel D). It is 
clear that similar residuals to those seen in the XTE J 1650-500 data 
at ~0.6-0.8keV are present, although they do not appear to be as 
prominent. This might indeed indicate the features arise at least 
in part due to calibration uncertainties. Alternatively, it might in- 
dicate they originate due to uncertainties in the absorption model 
used; were this to be the case the residual features observed in 
XTE J 1650-500 might be expected to be stronger than those in the 
Crab given the higher column density obtained for XTE J 1650-500. 
Therefore, we consider it likely that, at least to some extent, the soft 
residuals in XTE J 1650-500 are due to systematic uncertainties in 
either the XMM-Newton burst mode calibration or the absorption 
model. It is worth noting that residual features at soft energies have 
been seen in a number of timing/burst mode observations of XRBs 
which may also indicate calibration uncertaintie|f[ If we exclude 
the 0.7-1.0 keV energy range from our consideration, the simple 
reflection model (Model 1) gives a good fit, with xl = 1672/1601. 

In the interest of completeness, we also explore a number 
of physical explanations for the origin of these soft features, be- 
ginning here with interpretations that do not require additional 
absorption over that due to the ISM. However, we stress again 
our belief that systematic uncertainties are contributing towards 
the residuals, so the additional interpretations presented in this 
section should be treated with caution. We attempt to model the 
residuals as a combination of blended Ne II and Ne III absorption 
lines, whi ch are not includ ed in TBABS but frequently arise due to 
the ISM jjuett et al.l l2006l), and additional reflected OVIII emis- 
sion over that expected from an accretion disc with solar abun- 
dances, simila r to that seen in the ultra-compact neutron star XRB 
4U 0614-1-091 jMadei etalJboiOh . In the case of 4U06I4-H091 the 
strong oxygen emission was explained as arising due to accretion 
from an Oxygen rich white dwarf. Such a scenario is not plausi- 
ble for XTE JI650-500 so another source of oxygen rich mate- 
rial is required here. REFBHB is calculated assumin g solar abun- 
dances (as given in iMorrison &McCammoilll983l) so the extra 
OVIII emission is modelled with the addition of a narrow Gaus- 
sian at 0.654 keV to REFBHB before the relativistic effects were 
applied, while the Neon absorption is modelled simply with two 
narrow Gaussian absor ption lines, requir ed to contribute equally for 
simplicity (the work of lluett et alj|200d suggests this is not unrea- 
sonable). This combination provides a good representation of the 
data, with an improvement of 119 for 2 extra degree of freedom, 
and the parameters obtained are given in Table[3]as Model 4. 

On investigation, we find that either the combined Neon ab- 
sorption lines or the additional Oxygen emission can resolve the 
residuals individually, providing improvements of Ax^= 111 and 
1 14 (for one additional degree of freedom) respectively. For com- 
parison, the parameters obtained with these interpretations are also 
given in Table [3] as Model 2 and Model 3 respectively. How- 
ever, in the former case the Neon equivalent widths are EW = - 
55 ± 10 mA (for each line), much greater than those obtained using 
the 2001 Chandra MEG spectra of XTE J1650-500 (~ -2-10 mA; 
IMiller et al.l2004l ). In the latter, the disc inclination is strongly con- 
strained to be high (65 < i < 6 9°), inconsist e nt wit h that obtained 
with a similar interpretation in IMiller et alj l l2009l) . although we 



see |http://xmm2.esac.esa.int/docs/documents/C AL-TN-0083.pdf I 
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Table 3. Key Parameters obtained for the various models constructed using a reflection interpretation for the XMM-Newton spectrum of XTE J 1650-500 (see 
text); parameters marked with * have not been allowed to vary. Note that when included, the Nell and Nelll lines have been required to contribute equal 
amounts to the negative residuals at ~0.85keV. Column densities are given in 10^^ atom cm^^, abundances are quoted relative to the solar abundances, 
equivalent widths are given in eV, radii in gravitational radii, inclinations in degrees, temperatures in keV and hydrogen number densities are given in 10^*^ 
atom cm"^. For the model including the XSTAR giid, the ionisation parameter is given in erg cm s^^, and the outflow velocity in km s^^. Note that these 
results are obtained with unmodified XMM-Newton data; for results obtained after the application of the EPFAST tool see section !?. 1.71 
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Model 1 
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TBABS 
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'•^-0.1 
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, ^«4-0 04 

1 07^ 






, Qfi+o.oe 
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Ne II (abs) 
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400* 


400* 


400* 


400* 
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^^■^-0.2 


> 4.1 


7 3+1.1 
' "^-0.8 


> 4.3 




REFBHB 


kT 


0.200 ±0.001 
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1 1Q+0.04 
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•^■^-0.2 
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3 O+O-'^ 




O VIII (em) 


EW 






60+10 
0*^-30 


15+3 








1832/1663 


1721/1662 


1718/1662 


1713/1661 


1700/1660 



note that a different reflection code is used in this work, and the in- 
clination obtained here is consistent with t he lower lim i t of 50 ± 3° 
obtained from the optical light curve by lOrosz et al.l | |2004|) . We 
consider a combination of these features more plausible, and with 
this interpretation neither the inclination or the Neon equivalent 
widths are as tightly constrained, due to the degeneracy between 
the inclination and the emissivity index when modelling the line 
profile, and also the degeneracy between the individual interpre- 
tations. The Neon equivalent w i dths m ay be made consistent with 
those presented in iMiller et al ] 1 I2OO4I) by shifting more emphasis 
onto the Oxygen emission. However, we note that doing so tight- 
ens the requirement for a high inclination. It seems that in assign- 
ing a combination of these physical origins for the soft residuals 
this model can eith er be consistent with the inclinatio n obtained in 
lMilleretal.l l l2009l) or the Neon absorption presented in lMiller et al.l 
but it cannot be consistent with both. 

It is also possible that the soft residuals may indicate absorp- 
tion by ionised material associated with XTE J 1650-500 rather than 
the ISM, most likely due to an outflowing disc wind. To test this we 
generated a grid of ionised absorption models using XSTAR v2.2. 1, 
adopting solar abundances for all the elements and a turbulent ve- 
locity of 100 km s~i. The free parameters of this model are the 
column density, the redshift, which may in turn be used to infer 



the outflow velocity, and the ionisation parameter of the absorbing 
medium. Including this additional component in the basic disc re- 
flection interpretation constructed previously, the model also pro- 
vides a good representation of the data, with an improvement of 
Ax'^= 126 for 3 additional degrees of freedom. The key parame- 
ters for this interpretation are also given in Table[3] as Model 5. 

In this scenario, the residuals represent a complex blend of 
various narrow absorption features. The best fit value obtained for 
the redshift is positive, which suggests that if the material is as- 
sociated with XTE J1650-500 it is iiifalling, rather than outflow- 
ing, with a velocity of 6000 ± 4000 km s~i (although this is not 
particularly well constrained due to t he blend of multipl e features 
contributing). Given that the work of lOrosz et al.l ( |2004) strongly 
suggests XTE J1650-500 is a low mass X-ray binary (LMXB), 
and hence mass transfer should occur through Roche-lobe overflow 
rather than via stellar winds, the presence of such infalling material 
might be considered unusual, especially given the high velocity, so 
it is not clear this is the correct interpretation, although it remains 
an intriguing possibility. Such material could be the remnant of a 
failed outflow, which is now falling back onto the black hole. How- 
ever, even if the soft residuals are due to absorption by partially 
ionised matter, Fig. |4] shows that this material does not contribute 
any significant absorption features in the ~3-10keV energy range. 
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Figure 4. The model obtained treating tlie soft residuals as evidence for 
an ionised disc wind (see i]3.1.6l and Model 5 in Table [3) through use of 
an XSTAR absorption grid. The total model, corrected for Galactic absorp- 
tion, is shown in black, the red shows the contribution of the reflection com- 
ponent (including the thermal emission from the disc) and the blue shows 
the powerlaw continuum. Inset: the same model zoomed in on the 5-10 keV 
energy range. The ionised absorber does not contribute any strong absorp- 
tion features at these energies. 

In summary, a number of astrophysical origins for the soft 
residuals provide statistically acceptable fits to the data. However, 
given the physical issues each of these raise (requiring unreason- 
ably high neon absorption in the ISM and/or an unexplained over 
abundance of oxygen, or the presence of rapidly infalling material) 
and the comparison with the burst mode data from the Crab nebula, 
we strongly suggest these features are systematic rather than astro- 
physical. The roughly consistent inner radius obtained with each 
interpretation demonstrates that they are not critical to our determi- 
nation of this key parameter. 

3.1.6 Ionised Absorption in the Iron K Band 

We now focus on the 5-10 keV range in order to determine whether 
absorption from highly ionised species could modify the obtained 
line profile. The only possible absorption feature in the data is at 
~7.1keV. Assuming an identification with hydrogenic iron, this 
material would be outflowing at 6000 ±2000 km s~^, much higher 
than the typical velocities seen for XRB disc winds, which are 



less than ~1000 km s ^ (see e.^. 


Lee et al.l2002l,IUeda et alj2004 


iMiller et al.ll2006l,lNeilsen & Lee 


200S ), but not as high as the rel- 


ativistic wind described in Done & Gierliiiskil (|200S|) (see below). 



However, the detection is fairly tentative (including an unresolved 
Gaussian absorption line gives an improvement of Ax^~ 8 for 2 
extra degrees of freedom) and if real this is a weak feature (EW ~ - 
lOeV); its inclusion in any of the presented models does not alter 
the inner radius obtained. If this absorption line is real, and ionised 
absorption is the correct interpretation for the soft residuals, the sig- 
nificant difference in their velocities suggests these features cannot 
be associated with the same absorber. 

Based on the BeppoSAX observation of XTE J 1650-500 taken 
about a d ay before the XMM-Newton observation analysed here, 
iDone & G ierlifiski (2006) argue that a strong, highly blueshifted 
absorption line at ~7.6keV may be present due to ionised iron in 
a relativistically outflowing disc wind, and that failing to account 
for this line can lead to an artificially low inner disc radius. The 
line can instead be consistent with a truncated disc of inner ra- 




Energy (keV) 



Figure 5. iron line profiles for the XMM-Newton (black) and Bep- 
poSAX (red) observations of XTE J 1650-500, obtained on 2001 September 
13th and 11th respectively. The two observations may be fit with similar 
continuum models, resulting in practically identical line profiles. 



dius ~10-20f?G, although reflection of the Comptonised contin- 
uum by this disc, producing a fairly broad iron line, is still required. 
We extracted the BeppoSAX MECS spectrum in question (covering 
2.5-10.0 keV), and phenomenologically modelled the 2.5-4.0 and 
8.0-10.0 keV data with the basic BHB components as described in 
i]3.I.II the continuum model obtained for this observation is very 
similar to that for the XMM-Newton data. A broad excess emission 
feature is again clearly apparent when the 4.0-8.0 keV data are in- 
cluded. In Fig. [5] we show both the BeppoSAX and XMM-Newton 
line profiles; the two are found to be practically identical. 

We then verified that the MECS data, covering 2.5-10.0 keV 
can indeed be consistent with a truncated disc when ionised ab- 
sorption from a relativistic outflow is included, making use again 
of REFBHB and the XSTAR grid described above. Throughout 
this analysis a lower limit of 45° was imposed on the in clina- 
tion, similar to the lower limit found by lOrosz et al] j2004h . and 
the same interstellar absorption as obtained with the full XMM- 
Newton data was applied. We obtained an inner disc radius of 
16l3 for the MECS data, with absorber properties Nh = 
2.2lg:^ X 10^^ atom cm"^, ^ = 240 ± 60 erg cm s"^ and Wout = 
0.138 ± 0.007 c. These values are similar to those presented in 
iDone & Gie rliiiski' ('2006S. 

On investigation, even though there is no strong visual hint 
of any strong absorption in the high energy data, we also found 
that the 2.5-10.0 keV XMM-Newton spectrum can be consistent 
with a similar interpretation (i?in= 17tg^ -Rg), with iVH= 10^5 x 
10^^ atom cm"^, ^ = 140l^!^(, erg cm s"^ and Vont = S.Stto x 
10~^ c. However, as shown in Fig.|6](top panel), this absorber pre- 
dicts a variety of additional absorption features below 2.5 keV. In- 
cluding again the 0.5-1.7 keV XMM-Newton data, it is clear these 
absorption features are not present in the data, with the model pro- 
viding a very poor fit over this energy range 

(xl = 15477/1665) 

as demonstrated by Fig. |6] (bottom panel). Finally we attempted to 
model the full XMM-Newton spectrum with this interpretation, re- 
sulting in the inner radius settling again at ~2 Rg and the contribu- 
tion of the absorber being minimised. Although the interpretation 
invoking a combination of a truncated disc and a relativistic outflow 
remains viable if only the 2.5-10.0 keV data are considered, the ad- 
ditional data below 2.5 keV provided by XMM-Newton excludes the 
presence of absorption strong enough to affect the observed width 
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Figure 6. Top panel: The model obtained for XTE J 1650-500 considering 
only the XMM-Newton data above 2.5 keV and including a strong relativis- 
tic outflow, similar to that proposed in lOone & Gierliriskil 120061) . in order 
to allow the line profile to be consistent with a truncated disc. This out- 
flow would also imprint strong absorption features below 2.5 keV. Bottom 
panel: Data/model ratio plot of the broad band XMM-Newton spectrum to 
the model shown in the top panel, with the 0.5-1.7 keV data included. The 
low energy absorption predicted with the inclusion of a strong, relativistic 
outflow is cleaily not present in the data. 



of the iron line. Therefore, we consider the inner radius of ^2 Rq 
consistently obtained to be robust to the effects of absorption. 



3.1.7 Instrumental Ejfects 

As stated previously, the residual features seen in the ~ 1 .7-2.5 keV 
range are most likely due to calibration uncertainties associated 
with the silicon K and gold M instrumental edges. These fea- 
tures are most prominent for sources with high countrates, such 
as XTE J 1650-500, and it has been suggested that this may be (at 
least partly) due to residual charge transfer inefficienc}^ (hereafter 
CTI). This is a process that occurs during the readout of CCD de- 
tectors. As charge is transferred from one pixel to another, some of 
the electrons may be lost into 'charge traps' randomly distributed 
throughout the silicon lattice. This results in the measured energy 
of the incident photon being underpredicted. More recently, evi- 
dence has been found that suggests the residuals may instead be 
caused by X-ray loading (hereafter XRL). This effect is caused 
by X-ray photons contaminating the offset map produced prior to 
each observation to determine the energy zero-point for each pixel, 
and hence occurs preferentially for bright/hard sources. Whatever 
the physical cause, the current corrections applied in the standard 
EPIC-pn data reduction procedure do not adequately describe the 
instrumental features for high count rates. 

The EPFAST tool marks a serious effort to improve the per- 
formance of EPIC-pn fast-readout modes (Timing and Burst) 



''' see|http://xmm2.esac.esa .int/docs/documents/CAL-TN-0083.pdf| 
* http://xmm2.esac.esa.int/docs/documents/CAL-TN-0050-l-0.ps 




Energy (keV) 

Figure 7. Iron line profiles for the XMM-Newton (black) and Bep- 
poSAX (red) observations of XTE J 1650-500, as in Fig. |5] but here the 
XMM-Newton data has been modified with the EPFAST tool. This modi- 
fication results in clearly discrepant line profiles. 



by accounting for these remaining instrumental effects. In our 
reduction and analysis, we have found that EPFAST produces 
spectra that are better described by standard spectral responses 
at low energies, specifically in the 1.5-2.5 keV range where it 
corrects for the apparent shift of the instrumental edges very 
well. However, EPFAST is likely unsuited to higher energies 
at present, because it applies an energy independent correction 
{AE/E = constant). Studies show that CTI depends on energy 
as {AE/E)cTi oc i?"", with a ~ 0.5, meaning that the rela- 
tive impact of CTI is 2 times greater at 2 keV than at 6-7 keV (see 
e.g. lOendreau et alJll993llTownslev et alJl2002llGrant et alj|2004 
IPrigozhin et al.ll2004 iNaikaiima et alj|2008l) . Indeed, the CTI cor- 
rection included in the standard reduction procedures is an energy 
dependent correction, but any remaining CTI will also be energy 
dependent. If instead the residuals are due to XRL, the coiTection 
applied by EPFAST is an even worse approximation. To zeroth or- 
der, XRL causes a constant offset in the energy scale, i.e. AE 
= constant, therefore in this case, again writing {AE / E)xr'l oc 
E^°', we have a ~ 1, hence the relative impact of XRL at 6keV, 
in comparison to at 2 keV, is even less than for CTI. As the level 
of correction is primarily calibrated at ~2keV, especially for burst 
mode, EPFAST must over-correct spectra at high energies. 

This is supported by the striking difference between the Bep- 
poSAX and EPFAST-modified XMM-Newton spectra of XTE 1 1650- 
500 in the iron K band (see Fig.|7]l. In contrast, as demonstrated by 
Fig. [5] the unmodified XMM-Newton spectrum is remarkably simi- 
lar to the BeppoSAX data. As the XMM-Newton and BeppoSAX ob- 
servations in question are only separated by a day, given the long 
ter m evolution of XTE J165 0-500 during this outburst presented 
bv lDone & Gierlihskil ( [2006h . we expect any astrophysical changes 
displayed by XTE J1650-500 between these two observations to 
be minimal. The disagreement between the line profiles in the 
EPFAST-modified EPIC-pn and the BeppoSAX spectrum, as well as 
the agreement between the profiles in the latter and the unmodified 
EPIC-pn spectrum, suggest that the unmodified data may represent 
a better approximation of the true energy scale in the ~6 keV en- 
ergy range for count rates as high as measured in XTE J 1650-500. 
In an attempt to qualify the comparisons presented in Figs. |5] and 
|7]in a simple manner, we return to the Gaussian line profile con- 
sidered initially in section 13.1.21 and applied this model to both 
the EPFAST-modified and BeppoSAX data. The results obtained for 
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Figure 8. Left panel: x^confidence contours for the line energies obtained with the EPFAST-modified XMM-Newton (black) and BeppoSAX (red) datasets 
when modelling the iron emission with a simple Gaussian line profile. Horizontal dased lines show the 90, 95 and 99 per cent confidence intervals for one 
parameter of interest. The line energies obtained with the two datasets are clearly not in agreement. Right panel: Comparison of the XMM-Newton EPFAST- 
modified (black) and BeppoSAX (red) line profiles, with the energy scale of the former shifted by ~0.55 keV, as suggested by the left panel. Including this 
energy shift returns the profiles obtained to excellent agreement. 



these data, as well as those from the unmodified XMM-Newton data 
obtained previously, are presented in Table |4] As before, the line 
energies were constrained to be in the range 6.4-6.97 keV. 

The line width obtained with the EPFAST-modified data is 
marginally smaller than that obtained with the unmodified data. 
Although this must be considered a systematic change, as it arises 
due to a redistribution of photon energies, the difference is small 
in comparison with the statistical uncertainty on each; the widths 
of the two profiles are fully consistent within their 90 per cent 
confidence limits, and both values are extremely broad. Therefore, 
we conclude that in this case instrumental effects do not appear to 
be causing a statistically significant broadening of the line profile. 
However, when considering the line energies obtained, only the un- 
modified XMM-Newton data is in agreement with BeppoSAX . The 
disagreement between the EPFAST-modified XMM-Newton and the 
BeppoSAX data is further emphasised in Fig.[8](fe/i' panel), which 
shows the line energy confidence contours for the two datasets. 
The main effect of EPFAST in this case is to shift the line profile 
to higher energies, without significantly modifying its width. To 
further demonstrate this, in Fig.[8](ng/7f panel) we again overlay 
the EPFAST-modified and BeppoSAX line profiles, but now with the 
former shifted down in energy by ~0.55 keV in order to realign 
the line centroids; the profiles are returned to excellent agreement. 
Indeed, modelling the line with a LAOR2 profile at 7.5 keV (again 
shifted by ^^0.55 keV in comparison to the unmodified case) gives 
practically identical results to those presented in Table [T] i?in= 
l.'iltoxn Rg, ?innor > 8.1 and i = 70 ± 1°. Although a shift 
of 0.55 keV represents an energy difference of ^^8 per cent and EP- 
FAST only formally applies a gain shift of ~3 per cent in this case, 
in appendix A we demonstrate that the use of EPFAST naturally 
leads to a much larger energy shift when modelling a broad feature 
at 6.4 keV, fully consistent with that obtained, as the application of 
this tool modifies the high energy continuum, which can be difficult 
to determine underneath such a broad line. 

Based on this comparison and the theoretical considerations 
of the relevant instrumental effects discussed earlier, we conclude 
that the unmodified data do provide the best approximation of the 
true energy scale, while the modified data probably yield the best 
description of the line width. Furthermore, the good agreement of 



Table 4. Comparison of the results obtained for the unmodified and 
EPFAST-modified XMM-Newton calibrations and the contemporaneous 
BeppoSAX observation, when modelling the iron emission with a simple 
Gaussian line profile. The line widths are consistent in all cases, but only 
the unmodified XMM-Newton calibration provides a statistical agreement 
with the BeppoSAX data for the line energy (see also Fig. [8). 



Dataset 


-^Gauss 


cr 




(keV) 


(keV) 


BeppoSAX ( 1 1th September) 


< 6.41 


1.09 ±0.1 


XMM-Newton (unmodified) 


< 6.50 


^ 17+0.24 
-^■-^ ' -0.21 


XMM-Newton (EPFAST-modified) 


> 6.81 


1.07 ±0.25 



the BeppoSAX and unmodified XMM-Newton spectra suggests that 
some process other than CTI or XRL is the primary cause of the 
instrumental residuals at ~2 keV, the overall effect of this process 
at high energies (~ 6 keV) is small. We note briefly that if either 
CTI or XRL is the dominant cause of the residuals at ~2 keV, as 
appears likely, these effects could also be contributing to the soft 
residuals observed below 1 keV, given that the relative strength of 
both of these processes is higher at soft energies. These features 
are of little consequence for the main aim of this paper, but we 
again urge that the possible physical origins presented for the soft 
residuals be treated with caution. 

In the interest of completeness, and for further comparison, 
we also present the results obtained with the basic reflection inter- 
pretation (Model 1) applied to the EPFAST-modified data, includ- 
ing the 1.7-2.5 keV energy range. An excellent fit is obtained, with 
xt = 1843/1840. The values obtained for various key parameters 
are quoted in table |5] and we also list again the values obtained 
with the unmodified data for ease of comparison. We only consider 
the basic reflection interpretation here given the likely systematic 
origin of the soft residuals. In any case, with the EPFAST-modified 
data it is possible to resolve these residuals merely as absorption by 
the Galactic ISM, although some unusual abundances are required, 
potentially suggesting that the energy scale below ~2keV may also 
be incorrect. Our discussion henceforth is limited to the parameters 
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Table 5. Comparison of the parameter values obtained for the basic re- 
flection interpretation with both the unmodified and the EPFAST-modified 
XMM-Newton data. As with Table [5] column densities are given in 
10^^ atom cm^^, abundances are quoted relative to the solar abundances, 
radii in gravitational radii, inclinations in degrees, temperatures in keV and 
hydrogen number densities are given in 10'^'^ atom cm^'^. 



Parameter 



Unmodified 



EPFAST-modified 



Ao 

r 

^in 

i 
1 

kT 

"H 



rj -17+0.08 

'•J^'-o.og 

J^-fD_Q Q2 

^■'-'"-0.03 

95+0-04 
^•^''-O.OT 

n 0+0.2 
^•^-0.3 

55 ±2 
0.200 ±0.001 



+0.06 
0.04 



6.7 ±0.1 
1.36; 
1.7 ±0.1 
0.52 



+0.08 
0.04 



1 99+001 
3± 1 

0.200 ±0.002 



1.22 
5.0 



+0.03 
0.04 

+0.2 



1.34 
5.0 



+0.07 



+0.6 



that dictate the relativistic blurring, and hence determine the line 
profile: the main focus of this work. 

First of all, we note the consistency of the values obtained for 
the inner radius; the radius obtained with the EPFAST-modified data 
is marginally larger, but the statistical agreement with the value ob- 
tained with the unmodified data is excellent. Given the similarly 
excellent agreement of the line widths when using the simple Gaus- 
sian model, this is not surprising. There are, however, two signifi- 
cant differences that require highlighting. The inclination is signif- 
icantly higher and the index of the emissivity profile significantly 
lower when EPFAST is applied. This is a direct consequence of the 
line profile being shifted to higher energies by EPFAST. As we are 
using a physically self-consistent reflection model, the intrinsic en- 
ergy of the iron emission is hardwired into the model, and only the 
blurring kernel (KDBLUR) is able to change the energy at which the 
emission is observed. In order to shift the line to higher energies, 
the blurring kernel must both enhance the blue wing (increase the 
inclination) and attempt to suppress some of the red wing (decrease 
the emissivity index) of the line profile. The best fit inclination ob- 
tained using the EPFAST tool is unphysically high, implying that we 
are observing the source practically through the plane of the disc, 
although the 90 per cent confidence limit does extend to slightly 
less extreme inclinations. Therefore, the incorrect energy scale of 
the EPFAST-modified data at high energies can have a significant 
impact on key physical parameters when using physical rather than 
phenomonological models. We stress again that these changes are 
not physical, and are purely a result of the line profile being incor- 
rectly shifted to higher energies when the EPFAST tool is applied. 

To summarise briefly, we have investigated in detail the ap- 
plication of the EPFAST tool, which was designed to correct the 
instrumental features observed in EPIC-pn data around ~2 keV at 
high count rates. The correction applied by EPFAST does not - at 
the time of writing - have the correct energy dependence for ei- 
ther CTI or XRL, the two likely physical causes of these features, 
and although this does not lead to any significant changes in the 
width of the observed line profile in this case - which is ultimately 
the main focus of this work - it does result in an incorrect energy 




Figure 9. confidence contours for the optical depth (t) and electron tem- 
perature (kTe ) obtained for the Comptonising medium with COMPTT for 
the XMM-Newton observation of XTE J 1650-500. The contours shown rep- 
resent the 90 (black), 95 (red) and 99 (green) per cent confidence contours 
for two parameters of interest. 



scale at high energies. As demonstrated through comparison with 
a contemporaneous BeppoSAX observation and our simulations in 
Appendix A, EPFAST incorrectly shifts the line profile to higher 
energies. This energy shift can have important implications for key 
parameters when using physically self-consistent reflection models 
in which the line energy is not a free parameter We will proceed 
henceforth with the use of the unmodified XMM-Newton data. 



3.1.8 Alternative Scenarios 

Until now we have been considering the case in which the line 
arises due to reflection close to the black hole and is broadened by 
the relativistic effects present in such a region. This model there- 
fore requires that the accretion disc, which in this scenario is the 
reflector, extends in (or close to) to the ISCO. However, as stated 
previously, there is some debate over whether this is the case at low 
accretion rates (specifically in the low/hard state). Having demon- 
strated that the line is not artificially broadened by instrumental ef- 
fects, we now investigate whether the line can be primarily broad- 
ened by some other physical mechanism. The main alternative is 
that broadening via electron scattering (Compton broadening) dom- 
inates. This process gives a line width per scattering, a, of: 



a 

E 



(2) 



where E is the energy of the line and Tc is the electron temper- 
ature jPozdnvakov et alJll983h . In addition to broadening discrete 
features, Compton scattering also amplifies the energy centroid by: 



E„ 



En, 



= 1 + 4 



+ 16 



+ ■ 



(3) 



where Sobs and Eint are the observed and intrinsic line energies 
respectively. Considering the likely geometry and processes occur- 
ring in BHBs, Compton broadening could potentially occur in three 
separate media: the disc itself, the corona, or some outflowing disc 
wind. We will consider these three cases individually. 

First we consider the scenario in which the breadth of the 
line is primarily due to Compton broadening in the corona. Mod- 
elling the line with a Gaussian profile, we obtained a line energy 
of '--^6.4keV and a width of cr '-^ 1.1 keV (see sections [3. 1.21 and 
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Figure 10. confidence contours for the inner radius for each model presented in Table|3] The horizontal dashed lines represent Ax^ equivalent to the 90, 
95 and 99 per cent confidence intervals. Small values for the inner radius are significantly favoured in each case. 



13. 1.7b . As there is no obvious narrow component to the line profile, 
invoking Compton broadening as the origin of its width requires 
that the optical depth of the corona r > 1, i.e. the vast majority of 
the line photons must experience at least one scattering. Here we 
consider the simplest case, in which r ~ 1 and the line is broad- 
ened by a single scattering process. Substituting the values obtained 
from the Gaussian profile into equation |2] we find that, in order for 
this to be the case, kTc — 6 keV. If we wish to invoke multiple scat- 
terings, i.e. increasing the optical depth, in order to obtain the same 
line width, To would have to be lower still. However, we remind the 
reader that Comptonisation of the line in the corona reduces its ob - 
served equivalent width by a factor of (see lPetrucci et alj200l[) . 
so the case with r ~ 1 is the most reasonable. 

A corona with t ^ I and fcTe ~ 6 keV would produce a curved 
continuum in the XMM-Newton bandpass. However, in section 
13.1.31 we demonstrated that the continuum is powerlaw-like, with 
minimal curvature. To further demonstrate that such a continuum 
is not observed. Fig. |9] shows the confidence contours for r and 
kTc when modelling the high energy spectrum with a combination 
of COMPTT and a Gaussian emission line. A combination of r ~ 1 
and kTa ~ 6 keV is clearly excluded. In addition, considering equa- 
tion|3] Comptonisation of the line in such a corona should shift the 
line centroid to ~6.75 keV (assuming an initial energy of 6.4 keV). 
However, the line energy is constrained to be below 6.5 keV. There- 
fore Compton broadening in the corona does not produce a self 
consistent interpretation for both the line profile and the observed 
continuum. It is worth noting that lHiemstra et all bOIlh present a 
similar consideration for the high energy spectrum of XTE 1652- 
453, and also conclude that the line profile in that source cannot be 
explained by Compton broadening in the corona. 

Second we consider electron scattering in the disc itself. As 
previously noted, this can be an important process in BHBs given 
their high disc temperatures. Indeed, we stress again that the re- 
flection model used here, REFBHB, has been calculated specifically 
for use with spectra from XRBs, and self-consistently includes this 
process at an appropriate level for such sources. In order to demon- 
strate that relativistic effects are preferred over Compton processes 
in the disc, in Fig.[TO]we show the x^confidence contours for the 
inner radius up to 100 Rg for each model presented in Table [3] In 
addition, we also show the case of the basic reflection model but 



with the '-^O.V-l.OkeV data excluded, referred to as Model 1*, to 
demonstrate that the presence of the soft residuals has little effect 
on the value obtained for the inner radius with this model. Increas- 
ing the inner radius has the effect of shifting the emphasis from 
a combination of Doppler and gravitational broadening to a com- 
bination of Doppler and Compton broadening. Clearly a small in- 
ner radius, and therefore the former scenario, is favoure d at greater 
than 9 9 per cent confidence in each case. Furthermore, iReis et aU 
( |2009|) demonstrate that solutions invoking highly ionised reflec- 
tion at large distances from BHBs are not physically consistent with 
the expectation of a standard thin disc. 



Finally, we also briefly consider the possibility that the line 
is not produced by a reflection process in the disc, but is instead 
produced via reprocessing of the X-rays from the central source 
in an outflow, then broadened by ele ctron scattering as it e scapes 
from the outflowing medium (see e.g. iTitarchuk et al]|2009l) . Such 
a model may statistically reproduce the line profiles observed in a 
variety of neutron st ar and black hole XRBs. However, as shown by 
ICacketteta"l] ( l2010h . when considered in detail there are a number 
of issues with this 'windline' interpretation. We refer the reader to 
that work for the full discussion, and merely summarize the main 
points here. First and foremost, in order to reproduce the line pro- 
files observed the outflows must be optically thick. However, broad 
iron lines are observed in XRBs across a broad range in inclination, 
hence if such lines are produced in a disc wind these outflows must 
cover a large solid angle. This is in contrast with both theoretical 
simulations of disc winds and observational evidence for absorp- 
tion by such outflows; see section|4]for a full consideration of the 
(lack of) connection between broad emission lines and disc winds. 



In addition, the mass outflow rates obtained in lTitarchuk et al.l 

l2009l for the BHB GX 339-4 and the Neutron star (NS) XRB SerX- 
1 are both comparable with the Eddington mass accretion rate. 
However, in the observations analysed in that work both sources 
are known to be accreting comfortably below the Eddington limit, 
so the mass outflow rate required by the windline interpretation ex- 
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Figure 11. Compaiison of the iron line profiles for the XMM-Newton 
(black) and the latter two BeppoSAX (red) observations of XTE J 1650- 
500 (21st September and 3rd October; top and bottom panels respectively). 
Again, the profiles appear to have similar breadth. A quantitative compari- 
son is given in Table [6] confirming the similarity. 

ceeds the mass inflow rate due to accretior0. Furthermore, broad 
iron lines are observed in NS XRBs across at least two orders 
of magnitude in luminosity, therefore the high mass outflow rate 
would need to be independent of the accretion rate. This require- 
ment is further supported by the lack of significant evolution in the 
line profile with accretion rate dis played by XTE J1650-500 (see 
section [3T9t and XTE J1752-223 jReis et al.l201ol) . Finally, broad 
lines and high frequency quasi-periodic oscill ations (HFQPOs) ar e 
observed simultaneously in NS XRBs (see e.g. lCackett et alj2009t) . 
For any viable HFQPO formation mechanism, their presence dic- 
tates that the inner disc must be directly observed, rather than ob- 
scured by an optically thick outflow. In addition, iHiemstra et al.l 
l|) argue that it is difficult to reconcile a Compton thick outflow 
with a Compton thin corona. Taking all this into consideration, we 
conclude that the windline interpretation is not a viable solution for 
the presence of broad iron lines in XRBs. 



3.1.9 Line Profile Evolution 

If the line is primarily broadened by Comptonisation in either the 
corona or in a truncated disc, which moves in as the outburst pro- 
gresses through the transition from the low/hard to the high/soft 
state, some evolution in the width of the line profile is expected as 
gravitational effects become more important, and the temperature 
of both the disc and the corona changes (typical electron temper- 
atures are fcTL ~ 50-100 keV in the l ow/hard state and fcTc ~ 
30-50 keV in the high/soft state; see e.g. lMalzac & Belmonl2009l 



^ Although the authors of that work make a claim to the contrary, they in- 
con'ectly compare the mass outflow rate to the inflow rate for the high/soft 
state, despite the wind parameters being obtained for low/hard state obser- 
vations, in which the accretion rate is much lower 



Table 6. Evolution of the inner radius for the BeppoSAX and XMM-Newton 
observations of XTE J 1650-500 during its 2001 outburst, obtained by mod- 
elling the iron K Hne line with a LAOR emission line profile. 



Observation Date 


Instrument 


Rg 


1 1th September 
13th September 
21st September 
3rd October 


BeppoSAX 
XMM-Newton 
BeppoSAX 
BeppoSAX 


. c:9-l-0.02 
J^-3^_0.03 

1.63 ±0.04 
1.31 ±0.01 

-I 0-1+0.05 
-^■■'-^-0.04 



iGierliriski et"ai]| 19991 [Sunvaev & Titarchuklll980l) . In order to de- 
termine whether any such evolution occurs we now also consider 
the line profiles obtained for the latter two BeppoSAX observations. 
We model the continuum in the same way detailed previously, with 
the same standard BHB components, requiring the photon index to 
be in the range 1.7-3.0. In Fig. [TT|we show a comparison of the 
XMM-Newton line profile to the profiles obtained in the second and 
third BeppoSAX observations (top and bottom panels respectively), 
similar to Fig. |5] although here the strength of the iron line rela- 
tive to the continuum has changed, so the ratios have been scaled 
to match the continuum and the peak of the emission. Again, the 
profiles appear to have a similar breadth and shape in both cases. 

Since in these cases the shape of the continuum is signifi- 
cantly different from observation to observation, we also quantify 
the comparison by modelling the iron lines from all four observa- 
tions with a LAOR profile (given that the additional parameters in- 
cluded in the LAOR2 model could not previously be constrained we 
have adopted the simpler version, which assumes a single power- 
law emissivity profile). In this case, we model all four observations 
simultaneously, requiring only that the inclination does not vary 
between them (although the line energies are again constrained to 
be between 6.4 and 6.97 keV). The evolution of the key parame- 
ter in determining the width of the line, Ri^, is shown in Table |6] 
Although there is some slight variation between observations, the 
inner radii obtained are all very similar, confirming the visual sug- 
gestion that there is minimal evolution of the line profile during the 
progression of the outburst. This lack of evolution strongly favours 
the interpretation in which the iron line arises through reflection 
from the inner regions of a disc with a constant and stable inner 
radius. We also note again that a similar la ck of evolution is dis- 
played by XTE J1752-223 l iReis et alj|2010h . so such behaviour is 
not restricted to this individual source. 

3.1.10 Black Hole Spin 

As we have demonstrated that the line profile is not modified by in- 
strumental effects, and that relativistic disc reflection is strongly 
favoured over scenarios dominated by Compton broadening, we 
now present a quantitative measurement of the spin of XTE J 1650- 
500, replacing KDBLUR with the KERRCONV convolution model 
l lBrenneman & Revnoldsll2006h . KERRCONV applies the same rel- 
ativistic effects as KDBLUR, assuming a broken powerlaw emissiv- 
ity profile (as in equation[T), but directly includes the dimensionless 
black hole spin as a parameter, defined as a * = cJ/GM^, where J 
is the angular momentum and AI the mass of the black hole, so that 
it must be in the range [-1,1] where a* = I is maximal prograde 
rotation. Its other parameters are the indices and the break radius of 
the emissivity profile, the inclination of the disc and its inner and 
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Figure 12. confidence contours for the spin parameter measurements presented in Table|7] demonstrating the spin values presented in Table|2]are the global 
minima across the allowed range of values for a* . Horizontal dashed lines represent Ax^ equivalent to the 90, 95 and 99 per cent confidence intervals. 
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Table 7. Black hole spin parameters obtained by replacing KDBLUR with 
KERRCONV for the five models given in Table [3] The spin measurement 
for the basic reflection model is calculated both including (Model 1) and ex- 
cluding (Model 1 *) the 0.7-1 .0 keV data, to assess the effect of the residuals 
on the spin measurement. 



+0.02 
-0.01 

+0.01 
-0.04 



-0.12 
-0.002 



+0.03 
-0.08 



+0.01 
-0.02 



outer radii. In this model the inner radius is given in terms of the 
inner-most stable circular orbit, and we obtain the spin of the black 
hole by assuming the disc extends all the way down to this radius. 
We again make the simplification of a single powerlaw emissivity 
profile by requiring the two indices to be identical and assigning 
the break radius some arbitrary value, and the outer radius of the 
disc is again set as the maximum value accepted by the model. The 
black hole spin, the index of the powerlaw emissivity profile and 
the disc inclination were allowed to vary. 

In each case, the models including KERRCONV provided sta- 
tisticaUy equivalent representations of the data as those using KD- 
BLUR, without significant differences in any of the common pa- 
rameters. The spin measurements obtained are given in Table [7] 
and the confidence contours are shown for each model in Fig. 
1121 We again show the spin measurement and contour for the ba- 
sic reflection model with the 0.7-1 .0 keV data excluded (Model 1 *) 
to investigate the effect of the soft residuals on the spin value ob- 
tained with this model. Clearly the presence of these features does 
not significantly modify the spin obtained, as the values obtained 
including and excluding this data are consistent within their statis- 
tical uncertainties. In all cases the black hole is rapidly rotating, but 
maximal spin is ruled out at the 90 per cent level. To summarize. 



the overall constraint placed using REFBHB is 0.84 ^ a* ^ 0.98, 
although we note that if, as expected, the soft residuals are not as- 
trophysical the lower limit on the spin increases as only Model 1 
and Model 1* then provide valid estimates. 

The spin estimates presented here are slightly higher than 
the value of a * = 0.7 9 ± 0.01 presented for the same dataset 
in iMiller et"al] ( |2009|) . On investigation, we find that this is 
primari ly due to the dif ferent reflection codes used in these 
works; Mil ler et alj ilOO'^ make use of the CDID reflection code 
( IBallantvne et al.ll200lh . while here we use REFBHB, an adapta- 
tion of the more recent REFLIONX code which has specifically been 
calculated for use with BHBs. REFLIONX (and hence REFBHB) in- 
cludes a far more comprehensive range of ionised species for the el- 
ements with high cosmic abundance. These two codes use different 
energy ranges for the incident continuum, which may have an ef- 
fect on the reflection spectrum produced. Upon replacing REFBHB 
with a combination of CDID and various multicolour disc models 
we also obtain lower values f or the spin of XTE J1650-500, con- 
sistent with those presented in iMiller et alj ||2009|). More recently, 
the reflection code XILLVER taarcia & Kallmanll2O10h has further 
extended the range of ionised species and atomic transitions con- 
sidered, but results obtained with the use of that reflection code are 
typically in good agreement with those obtained with REFLIONX 
when considering relativistically broadened scenarios. Since there 
is no version of XILLVER calculated specifically for use with BHBs 
currently available, we do not make use of that code here. 



3.2 MCG-6-30-15 

Having demonstrated that the line profile observed in XTE J 1650- 
500 has an astrophysical origin, and that the data significantly pre- 
fer the interpretation in which relativistic effects dominate, we now 
consider the case of the active galaxy MCG-6-30-15. 



3.2.1 Continuum Modelling 

We began by modelling the spectrum of MCG-6-30-15 with the 
standard AGN powerlaw-like Comptonised continuum. However, 
the 0.5-10.0 keV X-ray spectra of active galaxies are rarely feature- 
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Figure 13. Data/model ratio plots for MCG-6-30-15. Top panel: ratio plot 
to the continuum model fit to the 0.5^.0 and 8.0-10.0 keV data (see text); 
the combined profile of the broad and narrow iron emission components 
is clearly shown. Middle panel: ratio plot with no contribution from the 
LAOR2 line to demonstrate the profile of the broad line. Bottom panel: the 
ratio plot when both LAOR2 and Gaussian emission lines are included in 
the model; the iron emission is now clearly well modelled. The spectra have 
been rebinned for clarity. 



less. Type 1 AGN commonly display addit i onal emission at soft 
energies (see e.g. Gierliiiski & Don i |2004 ICrummvetaDbOOd 
^vliniutti et al. 2009J. In order to phenomenologically account for 
this soft excess, we also included a black body emission compo- 
nent. Both of these were modified by neutral absorption. As MCG- 
6-30-15 is an extragalactic source, we initially allowed for both 
Galactic and intrinsic absorption with two TBABS components, one 
local compone nt and the other a t the redshift of MCG-6-30-15 
{z = 0.007749: Ipisher et al.|[T99l) . The Galactic column used was 
3.92 X lO'^" atom c m~'^, based on recent 21 cm measurements 
jKalberla et"al]|2005h . while the intrinsic column was free to vary. 
However, in the course of the modelling we found that the intrinsic 
neutral component was not required, and so was removed. 

As has been widely reported in previous work (see 
e.g. [B allantyne et alJ bOoH ; lYoung et alj|2005l : iMiller et allboOSl ; 
IChian g & Fabian IQUS), the spectrum of MCG-6-30-15 also dis- 
plays clear signatures of ionised absorption. To account for these 
feature s, we include a model f or the ionised absorber based on 
that of lChiang&Fabiaij l l201lh . In brief, this is a three-zone ab- 
sorber which als o includes a treatment of the iron-L dust absorption 
jLee et al.l200ll) . The absorber consists of a highly ionised, rapidly 
outflowing absorption zone (log^ ~ 3.82, Uout ~ 1900 km s^'^), 
and two moderately ionised, mildly outflowing zones (log^ ~ 2.3 
and 1.3, Vout ~ 150 km s~^). These are individually included us- 
ing the XSTAR photoionisation code, while the iron-L dust absorp- 
tion is accounted for with the inclusion of an additional TBABS 
component. As with XTE J1650-500, the 4.0-8.0 keV data were 
not considered so that the continuum was determined with minimal 
contribution from the iron line region. This model provides a good 
representation of the remaining data, with x2 = 1101/934. 



3.2.2 Line Profile 

As with XTE J1650-500, we now include the 4.0-8.0 keV data. Fig. 
[T3](top panel) shows the data/model ratio of the full MCG-6-30-15 
0.5-10.0 keV spectrum to the continuum model determined in the 
previous section. Again, it is clear that there is a large, broad excess 
in the ratio spectrum in this energy range, very similar to that seen 
in XTE J 1650-500. We again include a LAOR2 line to model this 
feature as a relativistically broadened iron fluorescence emission 
line from the inner regions of the accretion disc. The rest-frame en- 
ergy of the emission line was again constrained to be between 6.4 
and 6.97 keV, the outer radius of the disc was taken to be 400 Rg, 
and the rest of the line parameters were free to vary. However, in 
this case, the LA0R2 line did not completely account for the shape 
of the emission feature, due to the additional presence of a nar- 
row core to the iron emission, which can clearly be seen in Fig. 
[T3l (top panel). This was modelled with a Gaussian emission line 
at 6.4 keV (rest frame); narrow components to the iron emission 
in AGN can arise due to reflection from more distant, cold mate- 
rial, e.g. the dusty torus, or due to re-emission from photoionised 
clouds in the broad line region. With the inclusion of these compo- 
nents, the model provides a good representation of the data, with 
x2 = 1817/1725. The key parameters obtained are given in Table [T] 
and the data/model ratio plot is shown in Fig.[T3](bottom panel). In 
order to demonstrate the profile of the broad iron emission. Fig. 1131 
(middle panel) also shows the data/model ratio with no contribution 
from the LAOR2 line. 

This simple analysis demonstrates the similarity of the broad 
excesses observed in both the stellar mass BHB XTE J1650-500 
and the AGN MCG-6-30-15. In both cases these excesses are con- 
sistent with a relativistically broadened iron emission line, origi- 
nating from the inner regions of the accretion disc around a rapidly 
rotating black hole. We will now proceed to explore this similarity 
in more detail. 



4 DISCUSSION AND COMPARISON 

We have analysed the 2001 XMM-Newton and BeppoSAX obser- 
vations of XTE J 1650-500. In agreement with previous works on 
these observations, we find that this source displays a broad ex- 
cess around ~7keV over the standard disc - plus-Comptonisation 
model for BHBs (see e.g. iMiller et al.ll2002l ; iMiniutti et alj|2004l; 
iMiller et"al]|2009l) . The profile of this excess is not particularly 
dependent on the choice of continuum model. Such excesses are 
frequently observed in other non-quiescent BHBs with high qual- 
ity data, and are usually interpreted as iron emission lines excited 
through illumination of the disc by the Comptonised continuum, 
then broadened and skewed by the relativistic effects present close 
to the compact object. We have confirmed that such disc reflec- 
tion models can successfully reproduce the observed data, both 
phenomenologically and through the use of th e more physical re- 
flection model REFBHB jRoss & Fabiai]|2007l) . Although the ba- 
sic reflection model provides a good representation of the broad- 
band shape of the spectrum, there appear to be residual features at 
~0.8 keV in the XMM-Newton data (see Fig. [T). Through a brief 
analysis of the 2002 XMM-Newton burst mode observations of the 
Crab nebula we find that qualitatively similar residuals are present. 
We therefore expect that, at least to some extent, these features are 
due to systematic uncertainties in either the XMM-Newton burst 
mode calibration or in the ISM absorption model used, although 
in the interest of completeness we also explore a number of pos- 
sible physical origins (see section [3TT3}. The exact origin of these 
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Figure 14. Data/model ratio (left panel) and Data-Model residual (right panel) plots for the 2001 XMM-Newton observations of XTE J1650-500 (black) and 
MCG-6-30-15 (red), highlighting the broad iron emission lines. Residuals are in units of (normalised) cts s~^ keV"^; those for MCG-6-30-15 should be 
multiplied by 10^^. The spectra have been rebinned for clarity. Aside from differences in the blue wings which imply different disc inclinations, the line 
profiles are quite similar. 



residual features is not particularly important to our broader discus- 
sion of tiie sources analysed in tlris work, as even if these residuals 
are evidence for absorption by ionised material, this material does 
not have any significant impact on the high energy spectrum around 
the iron K shell transitions (see Fig.|4l(. 

In addition, we have also briefly analysed the long 2001 
XMM-Newton observation of MCG-6-30-15. These observations 
of XTE J 1650-500 and MCG-6-30-15 were selected because 
they represent some of the highest quality X-ray datasets avail- 
able for XRBs and AGNs respectively, and previous work on 
these sources invoking relativistic disc refl ection has indicated that 
both are likely to be rapidly rotating (see iMiller et al.ll2009l and 
iBrenneman & Revnoldsi200a for XTE J 1650-500 and MCG-6-30- 
15 respectively). We have demonstrated that even after accounting 
for the ionised absorption present in MCG-6-30-15, it also appears 
to show a broad excess at ~6-7 keV. Although we limit ourselves 
to a phenomenological analysis of this source, this is sufficient for 
our needs, and we again show that this excess can also be mod- 
elled as a relativistically blurred iron emission line arising in the 
inner regions of the accretion disc due to reflection of the intrinsic 
Comptonised continuum. 

However, these simple demonstrations are not the main 
focus of this work; reflection interpretations very similar to 
those constructed here have previously been proposed for 
both of these sources, using both these observations and oth- 
ers, and have been shown to be successful (XTE J 1650- 



500: iMiller et al 


I2OO2I; iMiniutti et alj 12004"; 'Miller et al.l 


200s ; 


MCG-6-30-15: 


Fabian et al. 


2OO2I; Fabian & Vauehar 


2003: 


Ballantvne et alj l2003l: IVau 


jhan (& Fabiad 12004; Reeves et al. 


20061; Brenneman & Revnoldsl2006l; Miniutti et al.l2007l). In broad 



agreement with these previous works, our analysis indicates 
the black holes powering both sources are indeed likely to be 
rapidly rotating, although our quantitative spin measurements 
for XTE J 1 650-50 are typically higher than that obtained by 
IMiller et al] ( l2009h . As previously stated, in the case of MCG-6- 
30-15 partially covering absorption has also been shown to provide 
a good representation of the observed ene rgy spectrum and variabil - 
ity when allowed sufficient complexity l lMiller et al.ll2008l , l2009h . 
In this interpretation, the curvature at ~4keV is due to the blue 
wing of a wide absorption trough rather than relativistically broad- 



ened emission. Ilnoue & Matsumotq l l2003h and iMivakawa et al.l 
( l2009l) propose similar interpretations, in which a relativistically 
broadened iron line is not required, and the variability is dominated 
by changes in the absorbers. 

Even with the quality of data available for MCG-6-30-15 the 
disc reflection/light bending and partially covering/variable absorp- 
tion interpretations remain statistically indistinguishable. Instead 
the main focus of this work is to present a comparison between the 
apparent profiles of the the excesses observed at ~6-7 keV in these 
two sources, which we argue favours the reflection interpretation 
for this feature in AGN. We stress that this is not a proof, merely 
a series of logical arguments and comparisons which we find com- 
pelling, and present here to the reader. However, this is not a state- 
ment against complex absorption in AGNs. There is clear evidence 
for ionised absorption in the example of MCG-6-30-15 used here. 
We are merely arguing that this process is not the correct origin for 
the feature at ~6 keV interpreted here as a broad iron emission line. 

It is clear from Table [T] that the values for the key parame- 
ter in determining the red wing of the iron line, Rin, obtained for 
XTE JI650-500 and MCG-6-30-15 are similar. In Fig.[T4]we focus 
on the line profiles of the two sources. A visual inspection shows 
that, aside from minor differences in the blue wings due to the dif- 
ferent inclinations at which we are observing the inner discs of 
these sources, the two profiles are indeed very similar. Specifically, 
we draw the readers attention to the similarity in the curvature at 
~4-5 keV, which is interpreted here as the onset of the relativis- 
tically broadened iron emission line. This curvature plays a strong 
role in determining the inner radius of the accretion disc, so after 
a visual comparison it is not surprising that similar inner radii are 
obtained. Given the similarity of the observed profiles in these two 
sources, it is natural to assume that they arise from to the same 
physical process. Since the two sources considered are separated 
by many orders of magnitude in terms of their black hole mass, this 
process must be mass independent, and is hence likely to be atomic. 
Indeed, we have already stressed that such excesses are commonly 
observed in AGN and BHBs. In addition, fairl y broad excesses are 
also observed in neutron star XRBs (see e.^. [ Cackett et alj |20I(]L 
iReis et al.ll2009l , ID'Ai et alj2oTol iLin et al.ll2010t) . 
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4.1 Complex Absorption 

We have already stressed that both disc reflection and partially cov- 
ering absorption provide statistically acceptable representations of 
the available data for MCG-6-30-15 (in the absorption case the cur- 
vature at ^^4-5 keV is due to the blue wing of an absorption trough 
rather than excess emission). Indeed, both of these interpretations 
are atomic in nature, and should therefore be mass independent. 
Given their relative proximity. Galactic XRBs are the ideal testing 
ground to determine the origin of these excesses. They are bright 
enough that the absorption may be very well constrained, hence we 
now consider the possible origins of, and observational evidence 
for, absorption in XRBs, with particular reference to XTE J1650- 
500. It is well known that the ISM contains a combination of neu- 
tral and ionised material. However, the Galactic column is not suf- 
ficiently high for the partially ionised species to significantly alter 
the X-ray spectra of Galactic sources; this has b een demonstrated 
with detailed stu dies of a number of sightlines jjuett et al]|2004l ; 
ljuett et al]|2006h . In addition, it is well accepted that this absorp- 
tion should not vary significantly from epoch to epoch. Consid- 
ering the source itself, although the masses of the black hole and 
its companion sta r have yet to be well constrained, the work of 
lOrosz et al.l ( |2004|) strongly suggests XTE J1650-500 is a LMXB, 
and hence accretes via Roche Lobe overflow rather than through 
stellar winds, so the X-ray spectrum will not suffer from obscura- 
tion due to such outflows. Many of the other XRBs (neutron star or 
black hole) that show evidence for a broad excess emission com- 
ponent at ■~6-7keV are also LMXBs, and hence should also be 
uncomplicated by strong stellar winds from their companions. 

The only other major source of material that could signifi- 
cantly obscure the source are outflows from the disc itself. Ra- 
diatively driven disc winds are expected to become increasingly 
important as black ho les approach (and pos sibly exceed) the Ed- 
dington Limit (see e.g. iKing & Pound3l2003h . Hydrodynamic sim- 
ulations of such winds (on all scales, from AGN accretion discs 
to early proto-planetary discs) suggest these outflows should most 
severely modify the observed X-ray spectra at high inclination 
angles {i.e. s hallow angles wit h respect to the surface of th e 
disc), see e.g. JProga et all (l2000l);|P roga & Ka llmai] (l2004l2002h : 
ISchurch et all ( l2009t) : lOwen et alj 72010): Si m et al] j2010t) . etc. 



The exact angle at which these winds become important must de- 
pend on the mass and the Eddington fraction of the black hole, 
which will in turn play in important role in determining the ionisa- 
tion structure and the mass outflow rate of the winds. For a 10^ Mg 
AGN at an Eddington fraction of Le = O.S. ISchurch et alj l [2009h 
find that the wind should start to appreciably modify the observed 
X-ray spectra at inclinations greater than ~60° . We obtain a simi- 
lar inclination for XTE J1650-500 yet find no evidence for a sim- 
ilar level of absorption, although we note that the Eddington frac- 
tion of XTE J1650-500 in this ob servation is likely to b e comfort- 
ably lower than that considered in ISchurch et al] ( |2009|) : assuming 
Mbh~ 5 Mq gives Le ~ 0.1. 

Orbital parameters obtained for Galactic XRBs through op- 
tical studies of the companion stars indicate they are observe d 
at a large range of inclination angles (see e.g. lOrosj l2003h . 
Although we do not find any strong evidence here, ionised 
absorption associated with outflowing material is seen rela- 
tively frequently in LMXBs, with notably strong examples be- 

'2007h and 



ing GROJ1655- 40 jMiller et alj|2006 



GRS 1915+105 jKotani et alj I2OOOI: iLee et alj l20oi 



Diaz Trigo et al 



Ueda et aU 



I2OO9I : iNeilsen & Lee|[2009l) . As expected, this absorption is more 
prominent in high inclination sources, and/or those observed at high 



Eddington fractions. Even at relatively high inclinations, the ab- 
sorption usually manifests itself as weak, narrow line features best 
studied with grating spectrometers, although XRBs with the high- 
est inclinations frequently show 'dips' in intensity, during which 
the absorption increases a nd, in some cases, can severely mod- 
ify the spectrum (see e.g. ISidoli et all l2002l iBoirin et all l2005l 
iDiaz Trigo et al.ll2006h . 

The general picture is that when present, disc winds from 
XRBs are more ionised (at a given inclination angle) than those 
in AGN. This might not be surprising as the discs around XRBs are 
expected to be more ionised than their AGN counterparts, owing 
to their much higher temperature. Indeed, IProga & KallmarJ 1 I2OO2I) 
show that the winds in LMXBs will be thermally driven from the 
disc, and should in general be optically thin. Therefore, in the vast 
majority of cases, these outflows do not dominate the observed 
shape of the X-ray spectra of BHBs. The fact that the continuum 
of XTE J1650-500 presented here is well modelled with the stan- 
dard BHB components (thermal disc at soft energies and Comp- 
tonised emission at hard energies) leads us to conclude that the X- 
ray spectrum of XTE J 1650-500 has not been significantly modi- 
fied by any additional complex absorption in this observation. Fur- 
thermore, in i]3.1.6l we have demonstrated that the XMM-Newton 
data exclude the possibility that the iron line profile may have been 
m odified by a strong, relativ istic outflow similar to that proposed 
bv iDone & Gierliiiskjl ( I2OO6I) . We briefly note that high mass X- 
ray binaries (HMXBs) also frequently display fairly strong absorp- 
tion f eatures {e.g. Cygnus X-1; see lHanke e t al. 2009, Ha nke et al.l 
I2OIOI) , but the absorption in these sources is expected to originate 
due to the strong stellar winds of their companions from which the 
BH accretes. 

In contrast to the evidence for ionised outflows in LMXBs, 
although the profile of the excess emission evolves with inclina- 
tion, its presence does not. The growing picture then is that many 
observations of XRBs show evidence for broad excesses at around 
~6-7 keV, despite these sources being observed at a large range of 
inclinations. The vast majority of these sources do not display any 
evidence for significant modification by complex absorption above 
that expected due to the ISM. These excesses cannot then be asso- 
ciated with the presence of disc winds, and must be intrinsic to the 



4.2 Common Origin 

If these excesses are not associated with complex absorption, as 
we have argued is highly likely, they must arise due to some emis- 
sion process. Limiting ourselves to atomic processes, as we wish 
to simultaneously explain the excesses seen in XRBs and AGN, 
the energies at which they are observed strongly suggest an asso- 
ciation with iron emission. Previous works have provided detailed 
comparisons of the various processes which may broaden one or 
more iron emission line s into an emission fea ture as broad as those 
observed in XRBs {e.g. iHiemstra et al ] |20Ilh . and frequently con- 
clude that, assuming the emission features are truly as broad as 
they appear, it is highly likely they have primarily been broadened 
by relativistic effects, which in turn requires the intrinsic emission 
to arise close to the black hole. Indeed, we have considered the 
main alternative for XRBs, namely that Compton broadening dom- 
inates, in some detail for the case of XTE JI650-500, and have 
determined that the scenario in which relativistic broadening dom- 
inates is strongly preferred (see sections D.I.8| and |3.1.9t . We stress 
again that some contribution to the line width from Comptonisa- 
tion is unavoidable in, but this is self-consistently included in the 
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REFBHB reflection model primarily used here. Furthermore, as we 
are searching for a common origin for the features in XRBs and 
AGN, we note that historically several authors have considered the 
possibility that Compton broadening may be the dominant effect in 
producing the broad lines observed in AGN, and ha ve concluded 
that such a scenario is not viable for these sources (see | Fabian et"al] 



1 19951: iRevnolds & Wilmsl 20001 : iRuszkowski et alJbOOol) 

We must therefore consider whether the excesses may be ar- 
tificially broadened in some way. It has recently been claimed that 
instrumen tal effects, such as p i le-up, may contribute to the broad- 
ening (see lYamada et al Jl 20091 : [Pone & D iaz Trigo"2010':'Ng et al.l 
I2OIOI) . However, detailed simulations by|Miller et al. (2010a) have 
shown statistically that pile-up most likely has the opposite effect, 
and should cause a narrowing of the features. Furthermore, we note 
that broad features similar to that observed here have been detected 
in XRBs and AGNs over a variety of luminosity states, and with 
a variety of different detectors, including gas based spectrometers 
which do not suffer at all from pile-u p. For a rece nt review on 
detections of broad emission lines see iMillej l l2007h . In the spe- 
cific case presented here, we can be confident that the spectrum 
of XTE J1650-500 does not suffer from pile-up effects, as the ob- 
servation was taken with XMM-Newton in burst mode, which is 
specifically designed to observe extremely bright sources. This has 
been confirmed through comparison of the line profiles obtained 
in this observation and the BeppoSAX observation roughly a days 
prior: Fig. [5] shows the two are practically identical, and the Bep- 
poSAX MEGS detectors are gas based and thus do not suffer from 
pile-up. In addition, this strongly suggests that other instrumen- 
tal effects, e.g. GTI, are not significantly modifying the profile of 
the line, as the profile appears to be independent of detector type. 
Therefore, we can conclude that the width of the line has a true 
astrophysical origin. 

All the evidence points towards the excesses in XRBs be- 
ing due to relativistically-broadened iron emission. We stress again 
that, even though complex absorption may initially appear to be a 
plausible explanation for these features in AGN, it is not a viable 
solution for XRBs. Indeed, reflection of the Comptonised emission 
by the accretion disc is a natural and unavoidable consequence of 
the presence of both physical components to the accretion flow, as 
long as some fraction of this emission is directed towards the disc, 
as is expected to be the case. Provided the disc extends in to the re- 
gions of strong gravity immediately surrounding the black hole, a 
widely accepted geometry for black holes with high accretion rates, 
the reflected emission lines will be broadened and possibly skewed 
by the relativistic effects inherent in such a region. The high cosmic 
abundance and fluorescent yield of iron indicate that its primary 
emission features, which are known to be between approximately 
6.4 and 6.97 keV (rest frame), should be the most prominent. 

The presence of both an optically thick accretion disc and 
some Comptonising region which produces high energy photons 
are observationally confirmed for both XRBs and AGNs (see 
e.g. iRemillard & McChntocklllool iDone et al.ll2007l lElvisll2010l 
etc.), therefore disc reflection does not require any additional com- 
plex emission regions over those known to be present. In addi- 
tion, gravitational light bending is a natural and unavoidable con- 
sequence for emission originating in a region of strong gravity 
in the general relativistic model, so the proposed explanation for 
spectral variability and steep emissivity profiles does not require 
any additional physics over known gravitational effects. Therefore, 
we argue that both disc reflection and light bending are present in 
AGNs. For the case of MCG-6-30-15, further evidence in favour of 
the presence of disc reflection/light bending interpretation is pre- 



sented by iTumer et al.l | |2004|) . who undertake a detailed analysis 
of the data obtained with the reflection grating spectrometer aboard 
XMM-Newton, and find strong evidence for a variable continuum 
but very little evidence for variable absorption, as would have been 
expected if the spectral var iability was driven by changes in the ab- 
sorbers. We also note that lRossi et al. 1 1 I2OO5I) find that the relative 
variations of the iron line and the continuum in XTE J1650-500 
favour a lightbending interpretation. 

In summary, since the presence of both accretion discs and 
Comptonising regions in AGNs is known, disc reflection is indeed 
a logical consequence, as are relativistic effects when the disc ex- 
tends close to the black hole. This is the most self-consistent in- 
terpretation for the broad excesses observed at ~ 6keV. Given that 
this process is also observed in XRBs, for which interpretations 
invoking complex absorption troughs are not viable, we conclude 
that the fundamental emission components observed from accreting 
black holes are a direct thermal component from the accretion disc 
(which appears in the UV for AGN), a Comptonised component 
extending to high energies, and a reflection dominated component 
arising due to irradiation of the disc by the Comptonised emission. 
Any intervening absorption must necessarily act on all these intrin- 
sic emission components. Other interpretations proposed for AGN 
add unnecessary complexities to the theoretical framework of ac- 
cretion in strong gravity, in discord with Occam's razor 



5 CONCLUSIONS 

With a simple spectral analysis, we have demonstrated that the high 
quality spectra from the 2001 XMM-Newton observations of the 
stellar mass black hole binary XTE J 1650-500 and the active galaxy 
MCG-6-30-15 both exhibit similar broad excesses at ~4-7keV. 
These are well modelled with relativistically-broadened iron emis- 
sion lines. Such features are commonly observed in both low mass 
X-ray binaries and active galactic nuclei. However, in the case of 
active galaxies there is still a debate over whether this feature arises 
through fluorescent iron emission or through complex absorption 
processes. In this work, we highlight the similarities of the fea- 
tures, herein interpreted as broad iron lines, in these two sources. 
By stressing that the general accretion process is mass indepen- 
dent, as optically thick accretion discs and Comptonising regions 
are present in both classes of black hole, we argue that the simplest 
solution for these apparent excesses in AGN is the same as that 
for their lower mass analogues. Galactic X-ray binaries, where we 
hope to have convinced the reader that the presence of relativistic 
disc lines is confirmed. 
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APPENDIX A: GAIN FITS AND BROAD LINES 

Here we present simulations demonstrating the effect of applying 
a constant multiplicative gain shift on the results obtained mod- 
elling features in the iron K energy region. We simulate data using 
the EPIC-pn response files from a simple model consisting of an 
absorbed powerlaw (Nh = 7 x lO^^atom cm"^, T = 2.15) and 
a Gaussian emission line with an equivalent width of 350 eV at 
6.4 keV. The exposure time for each simulated dataset was 1000 s. 
We then attempt to model this data with the same components, but 
with a multiplicative 3 per cent gain shift applied to shift the re- 
sponse to lower energies (equivalent to shifting the data to higher 
energies by the same amount, the approach taken by EPFAST). 
These values were chosen to best represent the case presented in 
the main body of the paper. This process is repeated 1000 times, 
for both narrow (a = 0.0 keV) and broad (cr = 1.1 keV) line pro- 
files (line widths were fixed throughout the modelling process, but 
continuum parameters could vary). In each case, the obtained line 
energy was obtained, and compared with 6.4 keV. The results are 
displayed in Fig. lAll For narrow line profiles (left panel), the line 
energies are tightly clustered around ~6.6 keV (^--^3 per cent differ- 
ence), as expected given the applied gain shift. This is because the 
feature is strong and narrow, and hence changes in the underlying 
continuum do not significantly effect the energy obtained. 

The same is not true when the broad line profile is considered 
(Fig. lAll right panel). Even though a gain shift of 3 per cent is ap- 
plied, we see that the energy obtained is systematically shifted by 
a much larger amount, and the line centroid obtained is typically 
~7keV (~9 per cent shift). This is because, unlike for a narrow 
line, the energy obtained for a broad line does depend somewhat 
on the underlying continuum model. The application of a constant 
multiplicative gain shift falsely increases the high energy curvature 
present in the data relative to the original continuum model; data 
simulated with a simple powerlaw model display visible, concave 
curvature (turning up at high energies) when the response is shifted 
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Figure Al. Distributions of the energy centroids obtained for both narrow {a = 0.0 keV; left panel) and broad (cr = 1.1 keV; right panel) Gaussian line profiles 
from the gain fit simulations described in the text. The lines were initially simulated at 6.4 keV, and the energy centroids were measured again after the 
application of a 3 per cent constant multiplicative gain shift, as applied by EPFAST for the case of XTE J 1650-500. The line centroids for the broad profile are 
systematically shifted by much more than the gain shift applied. The energy shift obtained for the broad line with the EPFAST-modified data in this work is 
indicated with an an'ow, clearly consistent with the simulated distribution. 
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Figure A2. Left panel: Data simulated with an absorbed powerlaw model using the XMM-Newton responses. The same model has been fit to the data, but with 
a 3 per cent multiplicative gain fit applied. In addition to the expected differences at the major absorption and instrumental edges, the gain fit clearly results 
in additional, artificial curvature in the continuum at high energies. Right panel: A Comparison of the effect of applying EPFAST and the XSPEC 'gain fit' 
command (see text). The two are found to introduce identical curvature at high energies. 



to lower energies, see Fig. |A2] (left panel). This additional curva- 
ture conspires to further increase the energy centroid obtained for a 
broad line profile. 

Finally, before we can relate this back to the shift in the line 
centroid energy observed when applying EPFAST, we must first 
confirm that EPFAST and the 'gain fit' command in XSPEC used in 
the above simulations do indeed have the same effect at high ener- 
gies. To do this, we apply the basic reflection interpretation (model 
1) to the unmodified XMM-Newton data with the 3 per cent gain 
shift used in the above simulations. Only the overall normalisation 
of the model is free to vary, all other parameters are fixed at the 
values in Table [3] Examining the data/model ratio (Fig. IA21 right 
panel, top), the same high energy curvature as highlighted previ- 
ously can be seen. Next we apply this model (with no gain shift) 
to the EPFAST-modified data, again only allowing the overall nor- 
malisation to vary. Examining the data/model ratio here (Fig. lA2l 
right panel, bottom), we see identical curvature. It is clear from this 
comparison that the use of EPFAST has the same effect at high en- 
ergies as the XSPEC gain fit command, also introducing additional 
curvature into the high energy continuum. We therefore conclude 
that the shift of ~0.55keV in the line energy obtained with the 



EPFAST-modified data in the main body of this paper, highlighted 
again in Fig. |A3l has been strongly amplified by the incorrect influ- 
ence of EPFAST on the continuum, and is a natural consequence of 
its application in this instance. 
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Figure A3. The unmodified (top) and EPFAST-modified (bottom) XMM- 
Newton line profiles; line energies obtained when using a Gaussian profile 
are shown by the dashed lines. The additional curvature at high energies 
introduced by EPFAST causes the centroid of the broad line profile to shift 
by a larger amount than the tool formally applies. 
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